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Introduction 

An analysis for predicting the buckling load of a 
compression-loaded orthotropic plate with a centrally 
located cutout is presented in reference 1. In this ref- 
erence, the classical two-dimensional analysis for deter- 
mining the plate buckling load is converted to an ap- 
proximately equivalent one-dimensional analysis by ap- 
proximating the plate displacements with kinematically 
admissible series. The accuracy of the one-dimensional 
analysis is a function of the number of terms included 
in these series. The objectives of this report are to de- 
scribe the one-dimensional analysis for a particular dis- 
placement series and to describe a computer program, 
BUCKO, associated with this analysis. The computer 
program is both easy and inexpensive to use. 

The main body of this report presents an overview 
of the basic assumptions used in the analysis, the 
derivation and solution of the governing equations, and 
a description of computer program BUCKO. The details 
cf the analysis are presented in the appendixes along 
with a user’s guide and sample problems. Symbols 
appearing in the main body and appendixes are defined 
in appendix A. The person interested in simply using 
the computer program needs to read only the main body 
of the report and the appendix containing the user’s 
guide and sample problems. The remaining sections of 
this report present the details required to modify the 
computer program. 

Analysis 

The assumptions made in the one-dimensional anal- 
ysis are outlined as follows. The plate is assumed to 
be a balanced, symmetric laminate with uniform thick- 
ness. The unloaded edges of the plate are simply sup- 
ported and the loaded edges are simply supported or are 
clamped. The cutout is centrally located with two mu- 
tually orthogonal axes of symmetry that coincide with 
the longitudinal and transverse axes of the plate. The 
geometry and the loading conditions used in the analy- 
sis are shown in figures 1 and 2, respectively. Loading 
is applied by either uniformly displacing or uniformly 
stressing two opposite edges of the plate. The plate has 
length L = 2c and width W = 26. The origin of the co- 
ordinate system is at the center of the plate. The x-axis 
is parallel to the loading direction and the y-axis is per- 
pendicular to the loading direction. The unloaded edges 
are free to expand in-plane in the y-direction. The plate 
is assumed to deform into a symmetric shape having 
mutually orthogonal axes of symmetry that also coin- 
cide with the plate axes. In-plane and out-of-plane dis- 
placements are represented by truncated kinematically 
admissible series which reflect the above assumptions. 

In reference 1 a one-dimensional formulation of the 
classical two-dimensional buckling analysis was derived 


following the Kantorovich method (ref. 2, pp. 304-327). 
The one-dimensional analysis consists of two parts: cal- 
culation of the in-plane stress distribution prior to buck- 
ling, hereafter referred to as the “prebuckling problem,” 
and calculation of the plate displacements and load at 
buckling, hereafter referred to as the “buckling prob- 
lem.” The results from the prebuckling problem are 
used with the principle of minimum potential energy 
and the Trefftz criterion (ref. 3, pp. 365-367) to obtain 
the stability equations for the buckling problem. 

Prebuckling 

As mentioned above, the in-plane displacements 
i the middle surface of the plate, u° and v°, are 
approximated by truncated kinematically admissible 
series. These series contain terms that are products of 
trigonometric functions in the y-direction and unknown 
generalized displacement functions in the x-direction. 
The series for the in-plane displacements are 
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where A is the loading parameter and N is the num- 
ber of trigonometric terms used. Substituting these 
series into the membrane potential energy functional 
and integrating over y reduces the potential energy to 
a one-dimensional form with the generalized in-plane 
displacements uo(x), ti 2 fc-i(z), an( * v 2 *-i( x) appear- 
ing as unknowns. The constant u 0 is a parameter 
which is selected to make the residual normal force on 
the unloaded edges of the plate equal to zero. The 
one-dimensional equilibrium equations are obtained by 
applying the principle of minimum potential energy. 
These equations constitute a system of simultaneous 
linear second-order ordinary differential equations with 
variable coefficients. The second-order system of equa- 
tions is solved for the generalized displacements and 
their derivatives, and these values are used to calculate 
the prebuckling stress distribution. The derivation of 
the one-dimensional prebuckling equations for N = 3 is 
presented in appendix B. The N = 3 analysis is used in 
computer program BUCKO. 

Buckling 

The results from the prebuckling problem are used 
to determine the buckling load of the plate. The 
prebuckling stress distribution is substituted into the 
functional corresponding to the second variation of the 
potential energy due to membrane and bending action 
of the plate. The out-of-plane deflection of the plate 




Figure 1. Geometry and coordinate system for a rectangular plate containing a centrally located cutout. 




(b) Uniform compressive edge displacement. 


Figure 2. Loading conditions. 
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middle surface w° is approximated by the truncated 
kinematically admissible series 

s 

w°[x,y) = ^u> 2 fc_i(z)cos [(2k- 1)^|] (2) 

k= 1 

where S is the number of trigonometric terms used. 
Substituting this series into the functional and integrat- 
ing over y produces a one-dimensional form in which 
the generalized out-of-plane displacements at buckling 
W 2 *-i(z) appear as unknowns. The one-dimensional 
stability equations are obtained using the Trefftz crite- 
rion. These equations constitute a system of simultane- 
ous linear fourth-order ordinary differential equations. 
The corresponding analysis for S' = 3 is given in ap- 
pendix C. 

Numerical Solution of the Differential Equations 

The prebuckling and buckling equations are solved 
numerically to determine the buckling load of the plate. 
The prebuckling analysis yields 2 N 4- 1 simultaneous 
second-order homogeneous differential equations and 
4iV + 2 boundary conditions. Some of the equations for 
the boundary conditions are nonhomogeneous. Com- 
puter program BUCKO uses a subroutine called PAS- 
VAR (ref. 4) to solve the prebucklirg differential equa- 
tions. This subroutine applies to both linear and non- 
linear systems of simultaneous first-order ordinary dif- 
ferential equations subject to two-point boundary con- 
ditions. To use PASVAR, the system of second-order 
differential equations is converted to the correspond- 
ing first-order system of equations. An example of this 
conversion is also presented in appendix B for the pre- 
buckling analysis associated with N = 3. 

The buckling analysis yields S simultaneous fourth- 
order ordinary differential equations. These differen- 
tial equations are linear and homogeneous and have 
4 S homogeneous boundary conditions. These homo- 
geneous equations and boundary conditions constitute 
an eigenvalue problem for the loading parameter A in 
equation (1). The smallest value of A for which a non- 
trivial solution exists is taken as the critical value of the 
loading parameter A cr . The corresponding critical load 
P cr is given by the value of the applied loading scaled 
by the critical value of the loading parameter. The sub- 
routine PASVAR cannot solve eigenvalue problems and 
is not used to solve the differential equations for the 
buckling problem. Central finite differences are applied 
to the buckling differential equations and the boundary 
conditions. This algebraic eigenvalue problem is solved 
for the critical buckling load of the plate. In the fol- 
lowing discussion, the algebraic eigenvalue problem for 
the buckling analysis is presented. Also, the modifica- 
tions required to use PASVAR in solving the prebuck- 
ling equations of the stress loading case are discussed. 


Algebraic eigenvalue problem for buckling . The 

differential equations of the one-dimensional analysis 
can be expressed as 

5 

[™2n-l(x)] 

n= 1 

- A//2 m -i,2n-l [^2n-lW] } = 0 

(m = 1,2, ... ,5) (3) 

where L 2m -i, 2 n-i (*) and ) are the differ- 

ential operators and W 2 n-i{x) are the generalized out- 
of-plane displacements at buckling. Using L T8 ( •) to rep- 
resent any L 2m _i,2n-i( ) operator, # rs (*) to represent 
any i/ 2 m-i,2n-*i (') operator, and iu(x) to represent any 
generalized buckling displacement, the form of these op- 
erators can be expressed as 

L r9 {w) = f\ (x)w{x) + [/ 2 (x)u/(x)]' + [/ 3 (x)u>(x)]" 

+ / 4 (x)w w (x) + [/ 5 (*K(*)r (4) 

and 

H r ,(w) = p 1 (x)w(x) + [p 2 (z)w'(l)]' 

+ [P3(*Mz)]' + p 4 (x)w'(x) (5) 

where /fc(x) and p*(x) are coefficients appearing in the 
differential operators and primes denote differentiation 
with respect to x. 

In deriving the finite difference expressions for these 
operators, the definition of the zth station derivative 
based on equally spaced half central differences is used, 
with the exception of the two terms in equation (5) 
multiplied by p 3 (x) and P 4 (x), The definition of the ith 
station derivative based on equally spaced full station 
central differences is used for these terms, since they 
are the only terms in equations (4) and (5) from which 
w x+l ^ 2 and w x ~ 1 / 2 unknowns would arise. Using this 
scheme for obtaining the finite difference expressions for 
higher derivatives and derivatives of composition func- 
tions, the discretized form of the differential operators 
given in equations (4) and (5) are 

[Lra(w))* = + to ’ -2 (/’- ‘/A 4 ) 

- 2 (/5 + / 5 -1 )/* 4 ] 

+ w’ [f\ - (/2 +1/2 + /2 -1/2 + 2/j + 2/j) /a 2 

+ (/^ +1 + 4/* +/*-»)/ A 4 ] 

+«•+> [(/• +i/ 2 +/r i +/i) / * 2 

- 2 (/; +1+ /s)/* 4 ] 

w+ 2 (ft 1 /* 4 ) 


(6) 



ltfr«(«0| l = + [?r 1/2 / a2 - (Po" 1 + P<) / 2A ] 

+*" , [p , .-(pi t,/2 +pr 1/2 )/ A ' 2 ] 

+ **+«[^ +l ' a /A a + (pJ +1 +Pi)H (7) 

where the superscripts refer to the finite difference sta- 
tion. Applying equations (6) and (7) to the operators 
appearing in equation (3) at each interior station and 
applying the basic difference expressions for the deriva- 
tives to the boundary conditions yields the algebraic 
eigenvalue problem 

\K]{w} = A[g]{u>} (8) 

The matrices [K] and [ g ] are square matrices of order 
S(M - 2), where M is the total number of difference 
stations. The vector {u/} contains the out-of-plane dis- 
placements at the interior difference stations. Because 
of the self-adjointness of the differential operators of 
equation (3), the matrices \K ] and [g] are symmetric. 
In addition, since the bending energy is positive defi- 
nite, the matrix [ K } is also positive definite. The small- 
est value of A for which equation ( 3 ) has a nontrivial 
solution is the critical value of the loading parameter. 

The finite difference and algebraic eigenvalue prob- 
lem formulations for the one-dimensional buckling anal- 
ysis corresponding to S = 3 are presented in ap- 
pendix D. 

Modifications for the stress loading case . The pre- 
buckling displacement field for the stress loading case 
can be determined only to within an arbitrary constant. 
This constant represents a rigid body displacement. In 
the present analysis, the Jacobian matrix for the pre- 
buckling differential equations (used in subroutine PAS- 
VAR) is singular because the displacement field is not 
unique. A procedure that corresponds to eliminating 
rigid body displacements from the displacement field is 
to reduce the number of prebuckling differential equa- 
tions and to rearrange the unknown generalized dis- 
placements. Appendix E describes these modifications 
to the prebuckling equations for the stress loading case 
for /V = 3. 

Description of Computer Program 

Computer program BUCKO is written in FOR- 
TRAN V for the Control Data Corporation (CDC) 
CYBER 170-series computers and is operational on 


the Network Operating System (NOS) level 1.4 at the 
NASA Langley Research Center. The computer pro- 
gram is composed of a main program, 10 subroutines, 
and 107 function subprograms, totaling approximately 
3600 lines of computer code. TV main program, pro- 
gram subroutines, and function subprograms are de- 
scribed in this section. 

The purposes of the main program are to read the 
input data file and to assemble the arrays used in 
the prebuckling and buckling problems. To minimize 
storage requirements, the arrays are stored in a master 
vector. The locations of the variables contained in the 
master vector are determined in the main program. 

The subroutine PREBUK is the major subroutine 
used in the prebuckling problem. This subroutine 
computes the approximate prebuckling stress distribu- 
tion. Subroutines JACOBI, JACOB2, FI, and F2 
are used in subroutine PREBUK. Subroutines JACOBI 
and JACOB2 determine the Jacobian matrices for the 
system of differential equations corresponding to the 
displacement loading case and to the stres : loading case, 
respectively. Subroutines FI and F2 calculate the co- 
efficients of the differential equations corresponding to 
the two loading cases. Function subprograms OMGl 
through OMG16 calculate the variable coefficients in 
the differential equations. The coefficients are functions 
of the cutout shape f(x). The cutout shape is deter- 
mined by the function subprogram R which is presently 
set up for elliptical and rectangular shapes and may 
be modified easily to include other shapes. Subroutine 
PREBUK calls subroutine PASVAR, which performs 
the actual numerical solution of the prebuckling equa- 
tions. 

The major subroutine used in the buckling solu- 
tion phase is BUCKL. This subroutine computes the 
finite difference coefficients, assembles the finite differ- 
ence equations, and computes the buckling load. Actual 
calculation of the eigenvalues is performed by subrou- 
tine SYMGEP, which is obtained from the NOS math 
ematics library. From the eigenvalues, the critical value 
of the applied loading is computed For the displace- 
ment loading case, the buckling load is obtained from 
the average normal stress acting at the loaded edges 
of the plate. For the stress loading case, the buckling 
displacement is taken as the average axial displacement 
of the loaded edges of the plate. Other output gener- 
ated by the program are the prebuckling stress distribu- 
tion, the buckling mode shape, the buckling coefficient, 
and average axial and transverse strains at buckling. 
A user’s guide and program sample problems are de- 
scribed in appendix F. 



Concluding Remarks 

This report presented an approximate analysis for 
predicting the buckling load of a rectangular orthotropic 
plate with a centrally located cutout. The analysis 
is applicable to plates that are compressed uniaxially 
by either uniform edge displacement or uniform edge 
Dormal stress. The boundary conditions considered 
are simply supported unloaded edges and clamped or 
simply supported loaded edges. A computer program, 
BUCKO, associated with this analysis was described. 


Buckling results generated by the computer program 
are values of the critical axial load, critical end short- 
ening, buckling coefficient, buckling mode shape, criti- 
cal average axial strain, and critical average transverse 
strain. 

Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
July 11, 1984 



Appendix A 
Symbols 

A 

All,A22, 

^12i^66 

b 

C 

! c\,c t , 

d 

Denl,Den2 

D 'i 

D\\,D22, 

D\2,D$g 

E \ 

/(*) 

fk(t) 

F 

Fb\,Fbz 


F x o, • • • 
F\ ^o, . . . 

■F tyl, . . 

Iff] 

[C] 


.Fts 

lFy5 

• > F x yh 


l G ]> G *i 

GbuGbi 

G) 

G x o. • • )G X 5 

GyO, . . • , G y5 
G*yl , • . . , G*y5 

Hbi>Hb2 

#;(■) 


cutout dimension defined in figure 1 

membrane stiffness coefficients of 
plate 

rectangular plate half-width 

rectangular plate half-length 

coefficients defined by equa- 
tions (C32) 

circular cutout diameter 

terms defined by equations (B59) 

finite difference coefficients defined by 
equations (D13), (D14), and (D15) 

bending stiffness coefficients of plate 


finite difference coefficients defined by 
equations (D4), (D5), and (D6) 

function representing cutout shape 
(see fig. 1) 

coefficients defined by equation (4) 

integrand defined by equation (C35) 

integrals defined by equations (C33) 

finite difference coefficients defined by 
equations (D7), (D8), and (D9) 

integrals defined by equations (C14), 
(CIS), and (C16) 


system geometric stiffness matrix 

system geometric stiffness matrix 
before application of the boundary 
conditions 

coefficients defined by equations (C7) 

integrals defined by equations (C33) 

finite difference coefficients defined by 
equations (D16), (D17), and (D18) 

integrals defined by equations (C8), 
(C9), and (CIO) 

integrals defined by equations (C33) 

differential operators defined by 
equations (C40) 




finite difference coefficients defined by 
equations (D19), (D20), and (D21) 


discretized differential operators 
defined by equation (7) 

H x 0i ■ ■ 

Hyo, ■ • 

H x yl, . 

.,Hx 5 
• > 

• • > Pxyb 

integrals defined by equations (Cll), 
(C12), and (C13) 

[K] 


system stiffness matrix 

[K\ 


system stiffness matrix before applica- 
tion of the boundary conditions 

L 


rectangular plate length 

MO 


differential operators defined by 
equations (C30) 

\L r ,(w 

)Y 

discretized differential operators 
defined by equation (6) 

M 


number of finite difference stations 

N 


number of terms in the prebuckling 
displacement series 

) Ny , N X y 

membrane stress resultants 

N° 


applied uniform stress loading 

N x0 ,.. 

iVyO, • • 

N X yl j • 

.,N X 5 

■ ,N y5 
. • , N X y 5 

stress resultants defined by equa- 
tions (B39), (B40), and (B41) 

Pk{x) 


coefficients defined by equation (5) 

Pb\,Pb\2.Pb2\ 

integrals defined by equations (C33) 

p ; 


finite difference coefficients defined by 
equations (DIO), (Dll), and (D12) 

p* 


applied axial force 

Py 


residual force on the unloaded edges 

PxCh • • 

P yOj • ■ 
Pxyl j • 

P yx 1 > • 

• » Pxb 

• i Py5 

• • > P ry 5 

• ■ > P yx 5 

integrals defined by equations (C17), 
(C18), (C19), and (C20) 

Qbi,Qbi 2 ,Qb 21 integrals defined by equations (C33) 

Qx 0) • • 

QyO, . • 
Qxyl i • 
Qyx 1) • 

• i Qi5 

• > Qyb 

• • j Qxyh 

• • > Qyx5 

integrals defined by equations (C21), 
(C22), (C23), and (C24) 

Rb1)Rb12, Rb2\ 

integrals defined by equations (C33) 

i?xo» • • 
RyOy • ♦ 
R x y \ , » 
Ryx \ » • 

• , R-o 
. , Ry& 

• * j Rxyb 
• • » Ryxb 

integrals defined by equations (C25), 
(C26), (C27), and (C28) 
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5 

number of terms in the buckling 

XI 

cutout dimension defined in figure 1 


displacement series 

x,y,z 

Cartesian coordinates 

Ub 

plate bending energy 

2/i,.--,2/u 

variables defined by equations (B42) 

u IS 

plate initial stress energy 


through (B49) 

U m , Urn 

membrane strain energy and mem- 

2/88 

a constant, see equation (E5) 


brane strain energy density 

2/o,.. -,2/i2 

variables defined by equations (El) 

u°, v° 

prebuckling displacements in the x- 
and y-direction, respectively 

7 

A 

displacement series parameter, 7r/2 b 
distance between equally spaced finite 

Wo, U2k-1, V2fc-1 

generalized displacements for the 
prebuckling problem 

6 

difference stations 
variational operator 

wo 

M 

displacement series parameter 
buckling displacement vector 


midplane membrane strains 
coefficients defined by equa- 

w° 

out-of-plane displacement 

tions (B58) 

w 2k-l 

generalized displacements for the 
buckling problem 

01) • • • >®23 

coefficients defined by equa- 
tions (B60) 

Wi,W3,W 5 

w 

w e 

values of the generalized displace- 
ments 

rectangular plate width 
external work 

A 

fil> ■ • • , ^16 

O' 

loading parameter 

integrals defined by equations (B6) 

differentiation of ( ) with respect to x 



Appendix B 

Prebuckling Analysis for N = 3 

The in-plane displacements prior to buckling are assumed to be 


u°(x<</) = A j«o(x) + £^U2*-i(x)cos((2fc- lby) j 
v°(x,y) = A jvoy + ^v 2 k-i(:r)sin[(2*- l)7y] j | 


BO 


where A is a loading parameter, 7 = tt/ 26 , and the functions uo(x), ttt(x), U3(x), «s(x), 17 (x), v 3 (x), and v 5 (x) are 
the generalized displacements to be determined. The parameter vo is a constant which is determined by requiring 
that the unloaded edges have zero resultant forces acting on them. Substituting these expressions into the linear 
strain-displacement relations from two-dimensional theory of elasticity yields 


e* = A juo(x) + ^U2 »-i(x)cos[( 2A — l)7y]| 

= A jv 0 + 7^(2* - l)v 2 *_i(x)cos[(2fc— l)7y)| 


Ixy 


=4ek*- 

U=i 


i(x) - (2k- l)7«2k-i'*)] sin |(2fc - l)7y) 


(B 2 ) 


where primes denote differentiation with respect to x. The membrane strain energy for the plate, which is assumed 
to deform symmetrically, is 

U m = f + ‘ f [An(e£) 2 + A 22 (e ° y ) 2 + 2 An(e° x e° y ) + ^6 6 (7* y ) 2 ] dy dx (B 3 ) 

J— c Jf(x) 

where f(x) is the curve shown in figure 1 and the Aij are the orthotropic membrane stiffness coefficients. Substituting 
the expressions for the strains into the membrane strain energy and performing some algebraic manipulations yields 

r+ c * 

U m = x 2 j f/ m (n 0 , Uo. u i. u 'i. u 3.«3. u 3 ,^,t> 5) ^) dx (B 4 ) 

where 


Um = [>lll (Uq ) 2 + ^22 V 2 0 + 2 Ai 2 UoI>o] ftl 

+ [■ 4 n(wi) 2 + + 2 Ai 27 u i v i] O5 

+ [i 4 n(u' 3 ) 2 + 94 3 27 2 t >3 + 64 i 27 u 3 u 3 ] fie 
+ [ 4 n(u' 5 ) 2 + 25i4 3 27 3 V5 + 10 j 4 i 27 «s w s] n i3 
+ 2 [/ 4 n uoui + Amv 0 vi + An ( 7 «oWi + vo«i)] 0 2 
+ 2 [AijUoUs + 3.4227^0^3 + An (37UqU 3 + vpu^)] 0 3 
-I- 2 [/llltloUi + 5 ^ 227 V 0 V 5 + ^12 (57U0V5 + U0U5)) Oxo 
+ 2 [An«l«3 + 3i4227 3 VjV 3 + Ai 27 (« 3 Wl d-Suivs)] O4 

+ 2 [i4iiU 3 tij + 15j4227 2t, 3V5 + -4l27(3«5 v 3 + 5U3V5)] flu 


Equation (B 5 ) continued on next page 
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fli(z) 

^(x) 

n 3 (*) 

n 4 (x) 

n 8 (*) 

n 6 (z) 

fl7(l) 

n 8 (x) 

n 9 (x) 

Dio(x) 

Hn(x) 

n»(x) 

ni 3 (*) 

n M (*) 

n l5 (x) 

n, 6 (x) 


+ 2 [AiiUjUg + bAvz^v iu 5 f Ai27( u 5 u i + 5 U 1 U 5 )] fii 2 

+ [^66 («i -7 «i) 2 1 0 8 

+ [^66( v 3 - 37U3) 2 ] fig 
+ [^66(^5 - 57U5) 2 ] n 16 
+ 2 [A66(vi - 7 tt i)( t 4 — 37 W 3 )] D 7 
+ 2 [^66(^3 _ $l u 3 )Ws - 57W5)] n 14 


[ dy-b- f(x) 
Jtt*) 

f b 1 

/ cos "iy dy = - \ 

JfW 1 

Z -6 

/ cos 373/ a;/ = - 
•'/(*) 

r b 




/(*) 

6 


/(*) 

f ' 

'/(*) 

rb 


-f, 


/(*) 

b 


/(*) 

6 


/(*) 

b 


-f 

Jf(x) 

-f 

Jf(x) 

-f 

Jf(x) 

-f 

•'/(*) 

-f 

■//(*) 

-f 

Jf l*) 

■f 

Jf(t) 


— j 

57 


- 57U5)] fils 

(B5) 


(B6a) 

- sin[7/(x)]} 

(B6b) 

{l + sin[37/(x)]} 

(B6c) 

= |sin[27/(x)] + ^sin[47/(x)] j 

(B6d) 

- /I 1 )] - ^ sin[27/(x)] 

(36e) 

“ /(^)) " sin[67/(x)] 

(B6f) 

= - jsin[27/(x)] - ^ sin[47/(x)] j 

(B6g) 

' /(*)1 + ^ sin[27/(x)] 

(B6h) 

- /(*)] + ^ sin[67/(x)] 

(B6i) 

1 -sin[57/(x)]} 

(B6j) 

= |sin[27/(x)] + isin[87/(x)] 

| (B6k) 

= “ sin[47/(x)] + ^ sin(67/(x); 

l} (B61) 

- /(*)] - 2^ sin[107/(x)] 

(B6m) 

= ~~ |sin[27/(x)] - i sin[87/(x)] j 

[ (B6n) 

= - 1 Pin[47/(x)] - isin[67/(x)] 

} (B60) 

" ^ + 267 8in I 10 ^/( jr )l 

(B6p) 
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When the loading is a uniform compressive edge displacement, no external work is done to the system. However, 
when the loading is a uniform compressive edge stress, that is, 


N x {±c, y) = -AJV* 


the external work is 


W e = -\N° X £ l«°( x ’V))lZV c dy 


(B7) 


(B8) 


The parameter N° is the applied stress shown in figure 2, and A is the loading parameter previously defined. 
Substituting the displacements given by equations (Bl) into equation (B8) for the external work yields 


<•+6 

w e = -x 2 n° J b («o(x) + Ui (x) cos 7 y 4- u 3 (x ) cos 37 y + u 5 (x) cos 57j/]*_^ dy 


The symmetry of the problem allows the external work to be expressed as 

rb 


By noting that 


(B9) 


W e = -2X 2 N° [ [uo(x) + tti(x)cos 7 j/ + u 3 (x)cos37j/ + « 5 (x)cos57y]* = ^ dy (BIO) 

Jo 


n,(c) = f dy = n 1 (— c) 

Jo 

Ojfc) = / cos 71/ dy = O 2 ( — c ) 

Jo 

n 3 (c) = / cos37j/ dy = fl 3 (-c) 

Jo 

ftio(c) = J cos 572/ dy - fi 10 (— c) j 


(Bll) 


the external work can be expressed as 

w e = -2X 2 n° ini (x)«o(x) + n 2 (x) Ul (x) + n 3 (x)u 3 (x) + n 10 (x)u 5 (x)]*z!^ 

The total potential energy of the system is 

V=--U m - W e 

which can be expressed as 

/ + C 

Urn (’io, « 0 > U 1 » U 1 > «3, « 3 , «5. Vi , v[ ,V 3 ,V 3 ,V S ,v’ & ) dx + 2X 2 N° (OiUfl + fytii + n 3 tl r + fli.0«5)*s 

(B14) 

The following ordinary differential equations and boundary conditions are obtained using the Euler-Lagr nge 
equations of variational calculus on equation (B14): 


(B12) 


(BI3) 


= +c 
c 
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tr 




Differential equations (on -c < x < -he) 


dUm d dU m 
duo dx \ du f 0 

dUm d fdUrn 
dui dx l du\ 


dUm d (dU^ 
du 3 dx \ du\ 


d (dU m \ 


dx y du' s ) 


= 0 






d . 

(dUm 


3t>i 

dx ' 

l dv[ 

= 0 






d . 

(dU m 


dv 3 

dx 1 

dv 3 

= 0 





dUm 

d 1 



dv s 

dx ' 

dv> 5 

= 0 



\ ° 


Boundary conditions (at x — ±c) 


2 thN? 6 uo = 0 


6u x = 0 


-?42n 3 AT* *u 3 =0 


— + 2n 10 ^ Uu 5 = o 


— 0 
6v 3 = 0 
6v 5 = 0 


Substituting the expression for U m given by equation (B5) into equations (B15) and 
Differential equations (—c < x < -he) 


(B16) gives 


^ [MllWo h ^i2t’o) ^1 -h (Au^i -h A\ 2lVi)tl2 


-h (4nu' 3 4 34^27^3) O3 4 (4n^5 4 54127^5) Oio] = 1 


4667 (K " 3 7«3) ^7 4 (v[ - 7«1 ) ^8 4 K - 57U5) n i5 J 


4 [(4hUq -h 4i 2 v 0 ) n 2 4 (Ann', 4 34i 2 7^3) H 4 


4 (A11U5 4 54i27 v 5) fli2 4 (4nu'j -h Ai 2 7 v i) ^5] = 0 


3A 6 e7 IK - 7 u i ) 0 7 4 ( v's - 37U3) ft 9 4 K - 57U5) n u ] 


4 ^ t(4nUo 4 4 i 2 v 0 ) n 3 4 (An *4 4* Ai27 v i) ^4 


4 CA11U3 4- 34 l2 7t'3) ^6 4 (Autig 4 5Ai27 v s) flu] = 0 


5A 6 g7 [(v[ - 71*1 ) f)i5 4 (t>3 - 371x3) n i4 4 K - 57u 5 ) Hi©] 


4 ^ [(AntiQ 4 Ai 2 v 0 ) fiio 4 (Anu'i 4 A127V1) fli 2 


4 (Anu'3 -r 34i 2 7u 3 ) fin 4 (A11U5 4 5 Ai 2 7 u 5 ) nj 3 ] = 0 


1 r- 
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7 [(AuUq + A 2 2 Vo) D2 + (A12U1 + A227V1) ft 5 

+ ( Ai 2 u ' 3 + 3^227^3) + (Ai2«5 + 5A227U5) n 12 ] 

- fa Mee [(v'l - 7«i) n 8 + (v' 3 - 37U3) n 7 + («5 - 57U 5 ) Di 5 ]} = 0 (B 21) 

^7 [(A12U0 + A22Vq^ O3 + (Ai2Uj + A227V1) fi 4 

+ (A12U3 + 3^227^3) n 6 + (A12U5 + 5A227W5) n n ] 

” ^ {^66 (K - 7 «i) D? + («3 - 37U3) n 9 + ( v ' 5 - 57U5) ni 4 ]} = 0 (B 22 ) 

57 [(Ai2«o + A22V0) Dio + (A12U1 + A227V1) fli2 

+ (A12U3 + 3A227V3) n„ + (A12W5 + 5 a 2 27 us) n 13 ] 

~ ^ {^66 (K - 7 «i) Dis + («>3 - 37 « 3 ) n 14 + (v' 5 - 57 u 5 ) n 16 ]} = 0 (B 23 ) 

Boundary conditions (x = ±c) 

{ [(-4n u o + ^12^0) Hi + (-AnUj -I- >li27Vi) n 2 

+ (A11U3 + 3Ai27V3)n 3 (AjiUg +5A127V5) n 10 + fIiiV°j 6u 0 } + =0 (B24) 

| [(AnUg + A12V0) D2 + {A\\u\ + Ai27t>i) Qs 

+ (A11U3 + 3Ai27V3)fl 4 + (A11U5 +5Ai27r 5 )ni2 +n 2 iV"j ^Ui| +C =0 (B25) 

{ [(Aiit*o + Ai2t> 0 ) fl 3 + (Aiju'j + Aj27v,)fl 4 

+ (Ai 1 u^ + 3A 1 27V3)n 6 + (A„u^+5A 1 27V5 )n„+n 3 N“] ^u 3 } +C = 0 (B26) 

| |(AiiUq + A12U0) Dio + (Anu', + A127U1) flj2 

+ (A11U3 + 3Ai27i» 3 )nn + (A,iu' 5 + 5A 12 7w 5 )ni3 + n 10 AT°J 6u 5 | +C =0 (B27) 

{[(u'l - 7«i) n 8 + (t>3 - 37U3) fi 7 + (t>5 - 57U5) Ois] < 5 t>i}-e = 0 (B 28 ) 

{K v i ~ 7 u i) d 7 + (i>3 - 37U3) n 9 + (V5 - 57U5) n 14 ] = 0 (B29) 

{(K -7“i)n 15 + («£ - 37 « 3 )n I4 +(v 5 - 57 u 5 )n 16 ] $v 5 }*® = 0 (B 30 ) 

The boundary conditions for the two loading cases investigated in the present study are given by 
Uniform compressive edge displacement 

u°{±c, y) = =f|a N xy (±c, y) rn 0 (B311 

Uniform compressive edge stress 

N x(±c, y) = -A N° N x v (±c, y) = 0 (B32) 
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Substituting the displacement series (eqs. (Bl)) into the displacement boundary conditions given in equa- 
tions (B31) and collecting like terms gives the following conditions: 

“ s(±c) = 0 1 {m 

ui(±c)=0 u 5 (±c)=0j 

Using the constkutive relations of two-dimensional orthotropic elasticity and the strain-displacement equations given 
by equations (B2) yields 


N x = A [(y4 n «o + A 12 vo) + {A n u\ + A 12 TfVi)cosTfy 

+ (A n u 3 + 3Ai27v 3 ) cos37y + {A u u' s + bA^v^) cos57y] (B34) 

N xy = AA 66 ((t>J - TfUi ) sin ~ty + (v 3 - Z^u 3 ) sin 3^y + (u 3 - 57115 ) sin 57y] (B35) 


Substituting these expressions into equations (B31) and (B32) yields the following conditions for the prescribed 
normal stress boundary condition in equations (B32), 

•An“o(±c) + A 12 V 0 = 

^nu'i( ±c ) +An'iM±c) = 0 
Anu' 3 (±c) + 3Ai2'yu 3 (±c) = 0 
Anu'^ic) + 5Ai27«s(±c) = 0 



and for the shear stress boundary conditions in equations (B31) and (B32), 

vj(±c) - 7 Ui(±c) = 0 
u 3 (±c) - 3 7U 3 (±c) = 0 - 
v' 5 (±c) - 57U 5 (±c) = 0 . 


(B37) 


Inspection of the conditions given by equations (B33), (B36), and (B37) reveals that these conditions are admissible 
since they satisfy the boundary conditions derived from the variational procedure. 

Upon solution of the system of ordinary differential equations for a specific loading case, the stress resultants can 
be obtained from the following equations: 


N x = X {N x 0 + N x 1 cos 7 y + N x3 cos 37y + N x5 cos 57y) 
N y = X (N y o + N y i cos 7 y + N y3 cos 37y + N y 5 cos 57y) 
N X y = A {N xy i sin 7y + N xy3 sin 37y + N xy 5 sin 57y) 

where 

N x o = AhUq + A12V0 ^ 

N xl = A u u\ + Anivi I 
N x3 = Anu 3 + 3A127V3 J 
Nxs = Aum's + 5A127V5 j 

Ny$ = Ai 2 Uq + A22V0 
Nyi = Ai 2 Uj + Aii'iv 1 
Ny 3 = AjaU 3 + 3^227^3 
JVys = A12U5 +5A227V5 , 


(B38) 


(B39) 


(B40) 
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N xy i = AeeM - 7^1) ' 
Nxy3 = AseK - 3 7“3) > 
W* v5 = AeeK - 57U5) , 


(B41) 


The differential equations given by equations (B17) through (B23) constitute a system of linear second-order 
homogeneous ordinary differential equations with variable coefficients, and these equations are solved numerically for 
the generalized displacements and their derivatives. To make use of existing computer subroutines, the second-order 
system is converted to an equivalent first-order system. The undetermined functions appearing in the displacement 
series are redefined as follows: 


«o(z) = yi 
U\(x) = 2/2 

U3(x) = 1/3 

“5(1) = V4 

In addition, the following functions are defined as 


t>i(x) = 2/5 
”3(2) = ye 
v 5 (z) = Vi 


(B42) 


y& — (Anj/J + Ai 2 t;o)Di -I- (Aiiy' 2 + Au'iys)^ + (^ 111/3 + 3^1271/6)^3 + (^ 112/4 + 5-4i27y7)Dio (B43) 

2/9 = CAnl/i + AizVo)tl 2 + (^ 112/2 + -^ 1271 / 5)05 + (A n y' 3 + 3^1272/6)^4 + {A\\y' A + 5Ai 2 72/7)ni2 (B44) 

2/10 — (-4nyJ + .4i2Vo)ft3 + (Anyj + A\ 2 iye)^A + (Any^ -I- 3Ai27ys)fl6 + (Any 4 + 5Ai27y7)fln (B45) 

2/ii = (•4 ii2/i + ^4i2Vo)Dio + (^nl/J + •4i272/5)fli2 + (■4n2/3 + 3 Ai272/6)Dij + (Any 4 + 5^1272/7)^13 (B46) 

2/12 = -466 [( 2/6 ~ 372/3) ^7 + ( 2/7 - 572/4) Di5 + ( 2/5 - 73/2) Dg] (B47) 

yi3 = -466 [(ys - TVi) B 7 + {y'e - 37y3) n 9 + (y 7 - 57y«) n M ] (B48) 

2/14 = -466 (( 1/5 ~ 72/2) fiis + ( 3/6 “ 37y3) On -I- (y7 - 57y4) flm] (B49) 

Substituting expressions (B43) through (B49) into the governing differential equations yields 

y'e = 0 (B50) 

y'g = - 7 yi 2 (B51) 

y '10 = -371/13 (B52) 

y '11 = -57yi4 (B53) 


7 [(-4, 2 yi + - 422 fo) H 2 + (Ai 2 y 2 + - 4227 y 5 ) D 5 

+ (-412V3 + 37-422V6) 0 4 + (Aial/i + 57-422y7) D 12 ] - y' 12 — 0 (B54) 

37 [(-4i2yi + A 22 V 0 ) D 3 + (-4i2yi + Ajiiys) fl 4 

+ (-4i2y3 + 3A 2 27ye) fie + (-4 13 y4 + 5A 22 7y7) Du] - y[ 3 = 0 (B55) 


& 


57 [(-4i2yi + -422^0) Dio + (A^yJ + Amye) fli2 

+ {Any 3 + 3A227V6) On + (Ai 2 yi + 5A227y7) r2i 3 ] - y' 14 = 0 (B56) 

The equivalent first-order system of differential equations is obtained by solving equations (B43) through (B46) 
simultaneously for yj, y' 2 , y' 3 , and y 4 ; solving equations (B47) through (B49) simultaneously for y' s , y' 6 , and y^; 
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and then substituting these values into equations (B54) through (B56) to obtain y' ]2 , y' 13 , and y' i4 . 
first-order system is given by 

.» _ , ?iys + ?2V9 + £31/10 + f4Vll 

5/1 " ^ Wo+ AnDenI 

A12 _ , £21/8 + f5l/9 + ffil/10 + £71/11 

»■ + 3^551 

/ ^12 , ^32/8 + ^62/9 + ^81/10 + ^11 

1/3 = a ^2/6 + 


y 4 = - 


An 

5Aj2 


7V7 + 


AnDenl 
£42/8 + ?7l/9 + ?9l/l0 + Sl02/ll 


An 

ys = iv-i + 

t/6 = 371/3 + 

y '7 - 571/4 + 


AnDenl 

^151/12 + ^162/13 + #172/14 

A66Den2 

^161/12 + ^182/13 + 0192/14 

A 6 6Den2 

^172/12 + #192/13 + 0202/14 


i4 6 6Den2 
yg = o 


y 9 = - 7 yi 2 

y'lo = “371/13 


where 


yii = -57V14 

y'i 2 = — — j “7 [^ 2^0 + 7 {ftsys + 3n 4 t / 6 + 5 r 2 i 2 i/ 7 )] + T^ 7 y 9 

An An 

yi 3 = v — “^7 [n 3 r 0 + 7 (n 4 y& + 3n 6 y6 + 50 n y7)] + -r-^yio 

All An 

fi 4 = " 1 - Y ^7 [nio^o + 7 (Diays + 3Hny6 + 5fii3t/7)] + -r-^TVn 

An An 


51 = fln07 + Oi2®8 + ^13^1 

52 = ^1302 — ^1109 ~ ^1008 
f3 = Oi 3 0 3 + Oi209 ~ Oio07 

£4 — -fiiotfi - o n 0 3 - n l2 fl2 

£5 = DiO^lO + Dij^ii -1- fii304 


£6 = fll0®12 + ^11013 + Oi305 
£7 = “01002 “ On05 - 0x204 

£8 = O 1O 014 “ Oi 2 013 + Oi306 ' 

£9 = “Oio03 “ On06 ~ Oi205 
£l0 — O 3 03 + 0405 + 0606 


Deni = 06021 + 03022 + 05023 “ O15 (0 2 Oio + O^n) ~ O10 (02O12 + 0lOio) 
Den2 = 02oOie + 0nOi5 + 0i9O; 4 


01 = 0 5 0 6 - Oj 

02 a O 3 O 4 — 0 a 0 6 

0 3 “ O2O4 — O3O5 

04 = o,o 6 -n§ 

0$ = O2O3 — O1O4 

06 * O1O5 - O a 

07 = O4O12 “ O5O11 

08 — O4O11 - OeOu 


09 — 0 3 0i2 ~ 0 2 0n 
0 io = 0 3 0u - OeOio 
0 n — 0 3 0io - OiOn 

012 = O 4 O, 0 - O3O12 

013 = O1O12 — O2O10 

014 = O2O12 - O5O10 

015 =* O 16 O 9 - 0?4 

016 = O14O15 - O7O16 


017 = O 7 O 14 — O 9 O 15 

018 = OeOie - 0i 5 

019 = O7O15 - 0 8 0u 

020 = OsOg — O 7 > 

021 = 0 6 0,3 - Oi! 

022 = 0 3 0i3 — 20ioOn 

023 = 040 } 3 — 20n0i2 


The resulting 

(B57a) 

(B57b) 

(B57c) 

(B57d) 

(B57e) 

(B57f) 

(B57g) 

(B57h) 

(B57i) 

(B57j) 

(B57k) 

(B571) 

(B57m) 

(B57n) 

(B58) 

(B59) 


(B60) 
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j The boundary conditions for the two loading cases are obtained from equations (B33), (B36), and (B37) and are 


j summarized as follows: 

Uniform compressive edge displacement 

y i(±c) = =f| 



(B61a) 


y 2 (±c) = 0 

(B61b) 


y3(±c) = 0 

(B61c) 


j y*(±c) = o 

(B61d) 


' 0i 5 (±c)yi 2 (±c) + tfi 6 (±c)yi 3 (±c) + 0i 7 (±c)yi 4 (±c) = 0 

(B61e) 


®i6(±c)yi 2 (±c) + 0 ig(±c)yi 3 (±c) + 0i 9 (±c)yi 4 (±c) = 0 

(B61f) 


j 0i7(±c)yi 2 (±c) + ^i 9 (±c)yi 3 (±c) + 0 2 o(±c)yi 4 (±c) = 0 

(B61g) 


Uniform compressive edge stress 



■\ fi(±c)ys(±c) +f 2 (±c)y 9 (±c) + f 3 (±c)yio(±c) + f 4 (±c)y u (±c) = -JV"Denl(±c) 

(B62a) 


• ? 2 (±c)y f (±c) + ft(±c)y 9 (±c) + f 6 (±c)y 10 (±c) + f 7 (±c)y n (±c) = 0 

(B62b) 


‘I ? 3 (±c)y tf (±c) + f 6 (±c)y 9 (±c) + ? 8 (±c)yi 0 (±c) + f 9 (±c)y n (±c) = 0 

(B62c) 


? 4 (±c)y 8 (±c) + f 7 (±c)y 9 (±c) + ? 9 (±c)yi 0 (±c) + fio(±c)y n (±c) = 0 

(B62d) 


0i5(±c)yj 2 (±c) + 0is(±c)yi 3 (±c) + 0i 7 (±c)yi 4 (±c) = 0 

(B62e) 


^i6(±c)yi 2 (±c) + tfi 8 (±c)yi 3 (±c) + 0i 9 (±c)y 14 (±c) = 0 

(B62f) 


®i 7 (±c)yi 2 (±c) + 0i 9 (±c)yi 3 (±c) + ® 2 o(±c)yi 4 (±c) = 0 

(B62g) 


For the uniform compressive stress loading, the differential equations and boundary conditions must be modified to 
make use of subroutine PASVAR. The modifications are given in appendix E. After numerical solution of the first- 
order system of equations, the stress resultants are given by equations (B38) and equations (B39) through (B41). 


In terms of the first-order variables, equations (B39) through (B4I) become 



AT fi ys + aye + ftyio + uvn 

Nx0 Deni 

(B63a) 


Ar _ fsys + isyo + &yio + fryn 
’ xl Deni 

(B63b) 


Ar _ &ys + f6ye + ftyio + ftyn 
J Deni 

(B63c) 


A7 uy& + fryd + &yio + fioyii 

^ l5 ~ Deni 

(B63d) 


AnA 22 “ ^12 , Ai 2 

NyO = v 0 + A N x0 

An A n 

(B63e) 


AT AnAa2-Aj 2 A 12 -- 

N y i= iys + An N * 1 

(B63f) 


1 n v3 = z AnA ”- A ' 2 ~m + ^N x3 

(B63g) 


I w,. = 5 4,,4 ?"' 4il Tvr + * , ’w.. 

J An An 

(B63h) 


*| 18 

1 


; 
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_ #15 3/12 + 0163/13 + 0172/14 

Den2 

N = - 1 — 12 + ^ 18S/13 + 0192/14 


tfxyl = 

^xy3 
^xy5 


Den2 

V _ 017J/12 + 0192/13 + 0203/14 

Jv iu5 — — 

Den2 


(B63i) 

(B63j) 

(B63k) 
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Appendix C 

Buckling Analysis for S = 3 

The prebuckling stress resultants used in this analysis are 


N x {x,y) = A j./V*o(x) + ^W i(2 *_i)(z)cos[( 2 *- l)iyj| 
N y {x,y) = A jiVyofa) + ^A r y (2fc-i)(x)cos[(2A:- l)-yyj | 
Vzy{x, y) = A|^ 7 V IS , (2fc _,)(x)sin[( 2 A:- lbj/]j 


(Cl) 


where A is a loading parameter, 7 = 7 r/ 26 , and the functions iV x0 (x), JV z i(x). A r X 3(x), ./V-^x), N y o(x), A r y i(x), 
N y 3(x), ./V^x), N xy i(x) 1 N xy z{x ), and N xy $(x) are determined from the prebuckling analysis given in appendix B 
(see eqs. (B 33 ) through (B 41 )). The series representation used for the out-of-plane buckling displacement is 


W°(x, y) = Y2 W 2j-1 (z) cos((2; - lj-yy] 

1 


(C 2 ) 


where the prescribed trigonometric functions satisfy both the kinematic and natural boundary conditions for a 
simply supported edge at y = ±b. The initial stress energy (contribution of the prebuckling stress distribution to 
the second variation of the potential energy) for the plate, which is assumed to deform symmetrically, is 


U,S = r [ b k* K) 2 + Ny (w° y f + 2N xy (»>•)] dy dx 

J~C Jf ( x ) 1 J 


(C 3 ) 


where f(x) is the curve shown in figure 1 and ( ) >x and ( ) >y denote partial differentiation with respect to the x- 
and y-coordinates, respectively. Substituting the series expressions for the stress resultants (eqs. (Cl)) and buckling 
displacement (eq. (C2)) into the initial stress energy (eq. (C 3 )) and integrating over the y-coordinate yields 


U is = A J {w} r [G]{tn} dx 


(C 4 ) 


where 


M - 


u>i{x)w[ (x)w f i(x) j ^(x^xju^x) j ^(xjiy^xju/^x) 


[Gn G 12 0 G14 G15 0 Gn G\% 0 

G22 0 G24 C25 0 G27 G28 0 

0 0 0 0 0 0 0 

G44 G45 0 G47 G48 0 
[ G ] — j G55 0 G57 Gsg 0 

0 0 0 0 

G77 G78 0 

G&g 0 

Symmetric 0J 


(C 5 ) 


(C6) 


Primes denote differentiation with respect to x. The nonzero expressions appearing in [GJ above are given by 
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i l '* 


Gil — NyoGyO + </Vyl< 7 j,l + Ny 3 Gy 3 4 Ny^Gy^ 

G \ 2 A^jryiGj-yj + N xy3 G X y 3 4 N X ybG X y§ 

G\A = NyoPyQ + Ny\Pyl + Ny 3 Py 3 *f Ny^Py^ 

G is = N xy iP yx i -f N xy 3 P yx3 4 - N xy 5 p yx5 
Gl 7 = NyoQyO + NylQyl 4 - N y 3 Q y 3 + N y ^Q y3 
^18 — ^ylQya-l 4* N xy 3 Q yx3 4- N xy 5 Q yx5 
G22 = N x 0 G x0 4- A^Gxi 4- N x 3 G x3 4- N x 5 G x5 
G24 — A^xyjPxyj 4- N xy 3 P xy3 4- N xy 3 P xy5 
G25 = N x 0 P x0 4- iVxiPxi 4 - N x 3 P x3 4* N x $Pzb 
G21 - N xyl Q xyl 4- N xy 3 Q xy3 4 - N xy 5 Q xy5 
G *28 — NxoQxO 4 A^xlQxl 4 N x sQx 3 4 N x3 Q x 5 
G44 = ^ y0 ify0 + JVyltfyl 4 JVy 3 tfy 3 + JVysffyS 
G 45 = N xyl H xyl 4 ” N xy 3 H X y 3 4* ATxysifxyg 

C47 = -NyO-RyO + A^yJ RyJ 4~ Ny 3 Ry 3 + Ny$Ry 3 

G 48 A^xyiPyxi 4“ ‘N I y 3 Pyx 3 4“ ^ xyb^yxh 
G55 — N x 0 H x0 4* iVxjPxl 4 N x 3 H x3 4- N x $H x s 

G 57 = -NjyjPxyl 4“ N X y 3 R X y 3 4- ATxy 5 P Z y5 

^58 = WxoPxO 4* iVxiPxl 4- N x3 R x3 4- N x3 R x b 
G 77 = NyoPyQ 4* A^yl^yl 4“ Ny 3 Fy 3 + Ny^Fy^ 

G 78 = -NxyjPxyl 4 A^xy3-^xy3 4 N xy $F X y$ 

G&b = WxqPxo 4 N x \F x i 4- N x 3 F x3 4- N X $F X 5 

where 

Gxo = / COS 2 7 !/ dy = 5 [<» - /(*)] - 7- sin[ 27 /(i)] 
y/(x) ^ 47 

Gl1 = /(x) C0s3 lydy= ^ 2 ~ sin M*)] { 2 + COs2 M*)]}) 

Gz 3 = J cos 2 7J/ cos 37 y dy = 4 (2 - sin [7/(1)) {2 + cos 2 [7/(1)] + 12 cos 4 [7/(1)]}) 

/* 6 _J 

Gl5 = 7/(x) C ° S2 iycos5iy dy = 420^ < 8 + 35sin [37/(1)] + 42 sin [57/(1)] + 15 sin [77/(1)]} 

Gy ° = Li*) 7 * S ‘ n2 7y dj/ = T {t 6- /(l) ) + ^ sin M*)) } 

/*6 

G »i = / 'Y 2 s^ 2 7V cos 71/ dy = 2 { 1 - sin 3 [7/(1)] } 

y/(x) »> 

/ 6 

( ^ 7 2 sin 2 7ycos37y dy = -^ (7 - sin 3 [7/(1)] {7 + 12 cos 2 [7/(1)]}) 


(C7a) 

(C7b) 

(C7c) 

(C7d) 

(C7e) 

(C7f) 

(C7g) 

(C7h) 

(C7i) 

(C7j) 

(C7k) 

(C71) 

(C7m) 

(C7n) 

(C7o) 

(C7p) 

(C7q) 

(C7r) 

(C7s) 

(C7t) 

(C7u) 


(C8a) 

(C8b) 

(C8c) 

(C8d) 


(C9a) 

(C9b) 

(C9c) 
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Gy5 ~ /(x) 72 Sin2 1V C0S &72/ dy = 4^0 {92 + 35 sin l 3 ^*)! - 42 sin [57/(*)I + 15 sin [7 T /(*)]} 


^iyl = 

^xy3 = 


^xy5 “ 


Z* 6 j 

- 7 cos 72/ sin 2 73/ dy = - {sin 3 [7/(2;)] - 1 } 

Jf(x) 6 

l*b J 

-J 7 cos 722 sin 7J/ sin 3 7 y dy = ± {sin [7/(2:)] - 1} - 1 {sin[57/(x)] - 1} 

f b 1 

- / 7 cos 72/ sin 7» sin 5 iy dy = — {4 + 7 sin [37/(1)] -3 sin [77/(1)]} 

Oft 


f b 1 1 

Hx0 = Jf(x ) COs2 37y = 2 ^ “ 12^ sin l 6 ^/( I )l 

y 6 j 

Fl1 = //(*) C0S7y COs2 37J/ dj/ = l40^ {72 " 5sin ! 7 T f /( z )] - 7»in(57/(*)] - 70sin[7/(x)]} 
ffl3 = L, , C0 ® 3 37y = i( sin [ 3 '7/( a: )H sin2 [ 3 7/(i)] - 3} - 2) 

f b 1 

^ lS = -//(X) C0s2 37yc0s57!/ dy = 2207^ 72 “ - 22sin[57/(x)] - 5sin[ll7/(x)]} 


Hs ° = /(x) 972 Sln2 37y ^ = ^ sinfSTT/ta:)] } 

f b g 

Hyl = //<«) 972sin237J/C0S ^ dy = ^{68 + 5sin[77/(x)] + 7sin(57/(x)] - 70sin[7/(x)]} 

^y3 = / 97 2 sin 2 372/ cos 373/ = — y{l sin 3 [37/(x)]} 

*'/(*) 

Hy * = /<.) 9l2ain2 ^ ycos5 ^ y dy = ^J{28- 55sin[7/(x)] + 22sin(57/(x)] - 5 sin[ll7/(x)]} 


H X y\ = 

H X y$ — 
HxyS = 


-J f{x) 3 7 sin 7y sin 372/ cos 372/ ^ ( 3 + J {5sin[77/(x)] - 7sin[57/(x)]}j 

~/( ) 3 'l ,sin23 'l r J' cos3 '7 d V= ^ {l + sin 3 [37/{x)]} 

y 6 _ 3 

■y 3 7 sin 3 7 j/ cos 372/ sin 572/ dy = —{12 - llsin[7/(x)) + sin(ll7/(x)]} 


F *° = J L °s 2 ^ y d V=2 [b ~ /{l) ) “ 20^ sm[107/(x)] 

/•b j 

Fxl = >/(,) C087l/CO82 ^ - 39^{ 300 - 188rin[7/(x)] - 11 sin[97/(x)] - 9 sin[ll 7 /(x)]} 

y 6 _j 

Fx3 = //(*) C0837yC082 57V ^ = l092^ {200 + l«28in[37/(x)] - 39sin[77/(x)] - 21 sin[13 7 /(x)] } 


(C9d) 

(ClOa) 

(ClOb) 

(ClOc) 

(Cl la) 
(Cl lb) 
(Cllc) 
(Clld) 

(C12a) 

(Cl2b) 

(Cl2c) 

(Cl2d) 

(Cl3a) 

(Cl3b) 

(Cl3c) 

(Cl 4a) 
(C14b) 
(Cl4c) 


APPENDIX C 


f b 1 

^ = 'iM COS35lV ^ = 157 “ 3sin (^/( :E )i +»in 3 [57/f*)]} 


(C14d) 


= / ;w 25f sin 1 5,, = S'! | 6 - + J_ Hl»lr/W]} (CIS.) 

" L) d « - S (294 “ 198 ™h/W] + Usin|9 7 /(x)] + »*,[U,/ W J} (CISb) 

fb 

^ = //(,) 25l2c0s3 ^ sin25 ^ d V = ^{-164- 182sin[3 7 /(x)] + 39sin[7 7 /(z)] 


FyO 

F 


'/(a) 

+ 2sin[137/(z)]} 

rb 


F vh - J 257 2 sin 2 57 y cos b^y dy = ^ {l - sin 3 [ 57 /(x)]} 


^iyl — 
^ry3 = 
-fiyS = 


rb 

-/^)57sm 7^57^71, dy = ^{30 - llsin[97/(z)] + 9sin[ll7/(x)]} 
/»6 r 

~J f{x) 57sin3 ^ sin5 ^ cos5 7y dy = A{20+ 13sin[77/(z)] - 7sin[137/(x)]} 

f b | 

- 57 sin 2 575 / cos 572 / dy = ^ {sm 3 [ 57 /(x)l - 1 } 


(C15c) 

(Cl5d) 

(Cl 6 a) 

(Cl 6 b) 

(Cl 6 c) 


P xO 

P« 1 

P*3 

PxS 


= cos 7 y cos 3 7 y dy = -i |sin[27/(x)] + 1 sin[4 7 /(x)]} {Cl7a) 

= JIM C ° S2 7J,COs37!/ dy = iJ^ ( X “ ^{3sin[57/(x)J + 10sin[37/(x)] + 15 sin[ 7 /(x)]} > ) (Cl7b) 

fb 1 ' 

1 W “■’»>» * - 55 ; (»- JfSsin^/WI + Tstal^/Wl + 70.toH/(x)])) (017c) 

f J 

— I cos 7 y cos 37y cos 57y dy 

Jim 

= 252^ {40 _ 63 s * n [7/( z )] - 21sin[37/(x)] - 9sin[77/(x)] - 7sin[97/(x)]} 


= / (i) 3 ^ sin ^8in37y dy = -i 7 {sin( 27 /(x)] - | sin[47/(x)] } 

/•fr ^ 

,yl = Jim ^ C087ysin ™ sin3 ™ dy = y (l + |{sin[57/(x)] - 5sin[ 7 /(x)]}) 

>V3 = Jim 37 ' sin7!/8in3 ^ C083 ^ i ( 3+ J{58in|77/(x)] - 7sin[57/(x)]}) 

rb J 

Pyi = I 37 3 sin 7 y sin 37y cos 57y dy 


P, 

P, 


>!M 

_ 1 


= 84 (-100 + 63sin[7/(x)] - 21 sin[37/(*)J - 9sm[77/(*)] + 7sin[97/(x)]} 


(Cl7d) 

(C18a) 

(Cl8b) 

(C18c) 

(Cl8d) 
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Pzyi - — f 37 cos 7 y 9in7j/ sin 37y dy = -^ ( 1 4- ^{sin[57/(z)] - 5sin[7/(z)]}') 

Jf(z) 5 \ 4 / 

^iy3 = ~ j f() 3 7cos7ysin 2 37y dy = ^17 + ^{5sin[77/(x)] + 7sin[57/(x)] - 70sin[7/(x)]}^ 

Pxy 5 = - 37cos7ysin37ysin57y dy 

JfM 

= -^{44-63sin[7/(z)] -21sin[37/(z)] +9sin[77/(z)] + 7sin[97/(z)]} 

Pyxi = - ^ 7 sin 2 7 yco 8 37y dy = ^7- ^{3sin[57/(x)] - 10sin[37/(z)] + 15sin[7/(x)]}^ 

Pyx3 = - fj^ 7 s *n IV sin 37y cos dy = ^3 + i{5sin[77/(x)] - 7sin[57/(z)]}^ 


yi5 


r 

= — I 7sin7ycos37ysin57 y dy 

Jf ( *) 


-1 


= 252(68 -63sin[7/(z)] + 21 sin[37/(z)] - 9sin(77/(z)] + 7sin[97/(x)]} 


Q*o = f cos 7y cos 57y dy = {3sin[47/(z)] + 2sin[67/(z)]} 

>'/(*) ‘47 

/ *> _j 

^ ^ cos 2 7ycos57y dy = ^^- {8 J - 35sin[37/(z)] + 42sin[57/(z)] + 15sin[77/(z)]} 

Qx 3= / cos 7ycos37y cos 57y dy 

•//(*) 

= 25^ 40 ~~ 63 sin[7/(z)] - 21 sin[37/(z)] - 9sin[77/(z)] - 7sin[97/(z)]} 

/*b ^ 

Q?5 - / cos 7y cos 2 57y dy = --{200 - 198sin[7/(z)] - llsin[97/(z)] - 9sin[ll7/(z)]} 
y/(i) oytry 

QyO = / 57 2 sin 7y sin 57y dy = ^{2sin(67/(z)j - 3sin[47/(z)]} 

J/(x) -64 

/*& 5-y 

Qy\ = / 57 2 sin7ycos7ysin57y dy = -57 {4 + 7sin[37/(z)] - 3sin[77/(z)]} 

Qy3 = / sin ^y cos 3^y sin 5^yy dy 

Jf(x) 

= ^{68 - 63sm[7/(x)] + 21 sin(37/(x)] - 9sin[77/(x)] + 7 sin[97/(x)]} 

Qys= / 57 3 sin7ysin57ycoG57ydy = ^{20- llsin[97/(x)]+9sin[ll7/(x)]l 

Jf(x) 


/*& 5 

Qx yi = - / 57 cos 7y sin 71/ sin 57 y dy = —{4 + 73in[37/(x)] - 3sin[77/(x)]} 

//(*) 84 


(Cl9a) 

(Cl9b) 

(Cl9c) 

(C20a) 

(C20b) 

(C20c) 

(C21a) 

(C21b) 

(C21c) 

(C21d) 

(C22a) 

(C22b) 

(C22c) 

(C22d) 

(C23a) 
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Qxytt = - I 5-7 cos 7V sin Zyy sin 57y dy 

<//(*) 


= 2^{44-63sin[7/(z)] -21sm[37/(i)] + 98in[7’7/(z)] + 7sin[9'y/(x)]} 


-5 

Qxy5 = ~ 57 cos 7y sin 2 572 / dy = — {196 - 198sin[7/(x)] + llsin[97/(x)] +9sin[ll7/(x)]} 

jf(x) 39b 


/’ fc —1 

Qyx i = - / 7 s >n 2 7y cos 572/ dy = + 35 sin[37/(x)] - 42sin[57/(x)] + 15sin[77/(x)l} 

•'/(x) 420 


Q y x 3 = ~ I 7 sin 72/ sin 372/ cos 572/ d2/ 
•'/(*) 


= ^2 {100- 63sin[7/(x)] + 21 sin[37/(x)] + 9sin[77/(x)] - 7sin[97/(x)]} 


f b -1 

Qyxs = ~ I 7 sin 72/ sin 572/ cos 57j/ dy = -{20- llsin(97/(x)] + 9sin[ll7/(x)]} 
d/(x) Jyb 


r 6 -1 

Rx o = I cos 372/ cos 572/ dy = — - {4 sin[27/(x)] + sin[87/(x)]} 
J}(x) 16 'y 

Rx i = t 
JfM 


cos 72/ cos 372/ cos 572/ dy 


= 252^ 40 “ 63 s» n ['7/(*) ] - 21sin[37/(x)] - 9sm[77/(x)] - 7sin[97/(x)]} 


i?i3 = / cos 2 372/ cos 572/ rfy 
•'/(*) 

1 


2207 ( 72 “ 55sin M*)l “ 22sin[5-y/(x)3 - 5sin[ll7,/(x)]} 

/*& _ | 

Rx 5 = cos 372/ cos 2 57y dy = - 

d/(x) 10927 


{200 + 182sin[37/(x)] + 39sin[77/(x)] + 21sin[137/(x)]} 


Ryo = I 157 2 sin37ysin57y dy = ^^{sin[87/(x)] - 4sin[27/(x)]} 

Jf(x) 16 

Ryi — / 157 2 cos 7y sin 37y sin 57y dy 

Jf(*) 

= ^-{44 ~ 63sin[7/(x)] - 21 sin[37/(x)] + 9sin[77/(x)] + 7sin[97/(x)]} 

Ry a •= f 157 2 sin 37y cos 37y sin 57y dy = -^{12 - llsin[7/(x)] + sin[ll7/(x)]} 

Jf(x) 44 

Ryi= [ 157 2 sin37ysin57ycos57y dy = —^{20 + 13sin(77/(x)l - 7sin[137/(x)l} 
•> f(x) 304 


R xy 1 = - / 578in7ycos37ysin57y dy 

<//(*) 


= {68 - 63 sin[ 7 /(x)] + 21 sin(37/(*)) - 9sm|77/(*)] + 7sin[97/(*)]} 


(C23b) 

(C23c) 

(C24a) 

(C24b) 

(C24c) 

(C25a) 

(C25b) 

(C25c) 

(C25d) 

(C26a) 

(C26b) 

(C26c) 

(C26d) 

(C27a) 
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f b -5 

R xy 3 = ~ 57 sin 371/ cos 37j/ sin 57J/ dy = ——{12 1! sin[7/(x)] + sin[l 17/(2)!} 

J}(x) 44 

R x y 5 = ~ f 57 cos 372/ sin 2 572/ dy = i“{164 + J82 sin[37/(x)] - 39sin[77/(x)] 

J f(x) lua * 

- 21 sin[137/(x)]} 

fb 

Ry X 1 = — / 37 sin 72/ sin 373/ cos 573/ dj/ 

•'/(*) 

= -^-{25 - 63sinb/(x)l + 21sin[37/(x)] + 9sin[77/9x)] - 7sin[97/(x)]} 

84 


Ry X 3 = - [ 37 sin 2 373/ cos 572/ dy = ;J -{28 - 55 sin[ 7 /(x)] + 22 s ; n[ 57 /(x)] - 5 sin[ll 7 /(x)]} 

7/(i) 

/* b 3 

i? VI 5 = ~ / 3 7 sin 3 7 V s * n 5 72 / cos &T 9 dy = — {20 + 13 sin[ 77 / (x)] - 7 sin[ 137 / (x)] } 

JfM 364 


'/(*) 

The strain energy of the plate due to the out-of-plane bending action is 

p+c rb 


U B = [ f \d u ( w° ) 2 + D22 (w° ) 2 + 2D 12 w° xx w° yy + 4 D 66 (w° „) 2 | dy dx 

J-c Jf(x) L J 

Substituting the buckling displacement series (eq. (C 2 )) into this expression and integrating over y yields 

r+c 


U B = j {m’} t [C]{w} dx 


where 



rc u 0 

C13 C14 0 C 16 

C17 0 C19 


<^22 

0 0 O25 0 

O 

s 

00 

O 



C33 C34 0 C36 

C37 0 C39 



C44 0 C4 6 

C47 0 C49 

[C] = 


c 55 0 

0 C 58 0 



<?66 

^67 0 c *69 




C77 0 C79 




c 88 0 


.Symmetric 


C99 

appearing in [C] are given by 


c n = 

D 22 Gbi 

^33 = D\\G X 0 

C55 = 4 Dee^yo 

Cl 3 — 

D\ 2 C B 2 

C34 = ^12^B12 

Csg = 4 D 6 qR v0 

Cu = 

D 22 P B \ 

C36 = £>11^x0 

C*66 = Z?ni/ x o 

Cl6 = 

DiiPbii 

C37 = D12QBI2 

Ce 7 = D\ 2 Hb \2 

Cn = 

D 22 Qbi 

C39 = £>llQ*0 

C *69 = D\\R X 0 

C 19 = 

D\ 2 Qb 21 

C44 = D22HBI 

C77 = D22FB1 

C22 = 

ADttGyQ 

C46 = X?12^B2 

C79 = D12FB2 

C 25 = 

*D W P y0 

C47 = D22RBI 

0 

so 

Q 

II 

oo 

C38 = 

4 DeeQ y o 

C49 = I?ia/?B21 

C99 = 


(C 27 b) 

(C 27 c) 

(C 28 a) 

(C 28 b) 

(C 28 c) 

(C 29 ) 

(C 30 ) 

(C 31 ) 

(C 32 ) 
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where 


G'bi — 


Gb7 

Hbi 

Hb7 

Fb i 
Fb i 

Pbi 

Fb 12 
Fb 21 
Rbi 
Rbi-2 
Rb 21 
Qb 1 
Qb 12 
Qb 21 


J ^'r 4 tos 2 ',i/ dy = y |[i> - /(x)] - ~sin[2'if/(i)]| 

= - J ^l 2 cos 2 nry dy = -y- | (6 - /(*)] - ^ sin(2-7/(z)j | 

= J ^ 8 17 4 cos 2 3^y dy = ^y- |[6 - /(x)J - ^ sin|67/(z)]| 

= - Jf 9"J 2 cos 2 3"yy dy - |[6 - /(z)] - sin[67/(z)] | 


-l 


625 7 4 cos 2 dy = 


2 

6257 




/(z) - — sin[10'y/(z)] 


} 


= - J 25^ 2 cos 2 5 ~iy dy = — y 7 - - /(z) - ^ sin(10-y/(z)! | 

— 9"y 3 

= / 97 4 cos 7y cos 37y dy = — - — {sm[47/(z)] + 2sin[27/(x)]} 

d/(a) 8 

= - / 97 2 cos 'vy COS 37y dy = — {sin[47 /(x) j + 2 sin[27 /(x) j } 

■'/(*) » 

A 6 i 

= - / 7 2 cos 7y cos 37y dy = ^{sinf47/(x)] + 2sin[27/(x)]} 

d/(z) 8 

r b — 2257 3 

= / 2257 4 cos 37y cos 57y dy = — — — {4 sin[27/(x)] + sin[87/(x] 

Jf{*) 16 

= - / 257 2 cos 37y cos 57y dy = ^{4sin[27/(x)] + sin[87/(x)]} 

Jf{x) 16 

f b O 97 

= - / 97 2 cos37ycos57y dy = — (4sin[27/(x)] + sin[87/(x)]} 
//W 16 

= f 257 4 


cos 7y cos 57y dy = 


-257 3 


{3sin(47/(x)] + 2sin[67/(x)]} 


/ 257 

— / 25-y 2 cos 7y cos 5-yy dy = — p{3sin[4'y/(z)] + 2sin[67/(z)]} 

Jf(x) 24 

/•* ~ 

- / 7 2 cos 7y cos 5^y dy = — {3sin[47/(z)] + 2 sin[6^/(z)] } 

.//(*) 24 


As in reference 1, the stability problem can be expressed by 

6(U IS + U B ) = 0 


(C33a) 

(C33b) 

(C33c) 

(C33d) 

(C33e) 

(C33f) 

(C33g) 

(C33h) 

(C33i) 

(C33j) 

(C33k) 

(C331) 

(C33m) 

(C33n) 

(C33o) 


(C34) 


This variational statement can be expressed as a functional containing the unknown functions in the buckling 
displacement and their derivatives, that is, 


/+« 

6 (Uis + U b ) = 6 J F (u>i, to,, w",w 3} w' 3 , w 3 , tu 5 , W 5 , 105 ) dx - 0 


(C35) 
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The following ordinary differential equations ?nd boundary conditions are obtained using the Euler- Lagrange 
equations of variational calculus on equation (C35): 

Differential equations (on —c < x < c) 


dF d / dF\ d?_ / 3 F \ 

dwi dx dx 2 V dw” ) 

dF d / dF \ <P_ ( dF \ _o 

dwz dx \dw 3 y dx 2 \ dw^ ) 

dF d ( dF _ \ <P_ ( dF\ 

dw$ dx + dx 2 \ditf 5 / 


Boundary conditions (at x = ±c) 


dF 

d dF 

dw[ 

dx dw'l 

dF 

d dF 

dw 3 

dx dw 3 

dF 

d dF 

dw' 5 

dx dw 5 


6w\ = 0 


SH- 

(5)**- 

(SH- 


Alternately, the system of ordinary differential equations can be written in operator form as 


L n ( ) L 13 ( -) i, 5 ( )l [//„(•) ff 13 ( ) ( tm(z) 

£3i(-) i 33 ( ) L 35 ( ) \ w 3 (x) > = A H 3 ,(•) tf 33 () i/ 35 ( ) < w 3 (i) 

Isi( ) Ls 3 (-) L 55 ( ) UsOOJ (•) H 53 ( ) ffss( ) [ «*(*) 


W) 

^i 3 () 

^15(‘) 

L* i() 
L 33 (-) 

i 3 5() 

^si(-) 

^55(0 


[Ci 3 (') + c 33 (-)T 

[C 34 () + C 36 ()T 
(G 37 () + G 39 ()"]" 
[C 16 (-) + C 36 ()"]" 
(C 46 () + C 66 ()"]" 
|G 67 () + C 69 ()T 
[C 19 () + C7 39 ()T 

[Cof*) + C69()"]" 
[C 79 (-) + C 99 ()T 


[C22(-)T + [Cn(') + £i 3 (T1 

[C25(-)T + [^h(-) + C 16 (-n 

[C2 8 ()T + (Ci7() + C 19 ()"] 

(C25(n , +[c ] 4()+c 34 (-n 
[c 55 (-)r + [C44 (*) + C46()”] 
[C58()T + (^7(-) + C49(-n 
[C28 (’) , ] / + [<^17 (*) + G 37 (-)"] 
[C 5 8(-)T + (C47(0 -I- <^67( )"] 
[C88(-)T + (^77(0 + <? 7 9(T] 


Hiz() 
»!«(•) 
Hz i() 
Hz 3 (-) 
Hzz(-) 
H*i() 
Hzz(-) 
Hz 5 (0 


[Gi2(0 + G22()T 
[G24 (') + G 2 5(-)T 
[G?2 7 (*) + C?28 ( ' ) , ] / 

(Gis(-) + G 25 (-)T 

(G45(0 + G 5 5(')']' 

[Gs 7 (-) + GseCOT 
[Gi8(0 +G28(-)T 
1^48(0 + G 88 (')T 
t<j r 78() + G88(')T 


[Gn(-) + Gi 2 (-)'] 
[Gu(-) + Gi 5 (-)'] 
(Gi 7 (-) + Gi8 (-) , J 
(Gm(-) + G 24O)'] 
(G44(-) + G 45 (')'] 

[g 47 (-) + g 48 (-)'] 

(Gi 7 (-) + G 2 7 (-) / ] 
IG 47 ( •) + G 57 (')'] 
[G 77 (-) + G 7 8 ( 0 '] 


A 
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Finite Difference Equations for S = 3 


For the case of S = 3, the differential equations of the buckling problem are 


3 

T: {^2*-1,2j- 1 ( w 2j-l(z)] - ^#2t-l,2j-l [w2j-\ (*))} = 0 (* = 1, 2, . . . , 5) 

J = l 


(Dl) 


Expressions for the differential operators are given by equations (C39) and (C40). Applying the difference expressions 
for the operators given by equations (6) and (7) to equations (C39) and (C40) yields the following finite difference 
equations (superscript i refers to the tth difference station): 


|L„ («,,)]♦' = E\w \~ 2 + E\w\~ l + + E' 4 w\ +l + E\w \ +2 

[£«(«*)]* = Eiw \~ 2 + E' 7 w' 3 -' + E' e w ' 3 + + E\ 0 w ' 3 + 2 

[LjsM* = E' n wl ~ 2 + E\ 2 w'f l + E\ 3 wi + E\ 4 w \ +1 + E\ s wi +2 
(I 3 iK)]‘ = F{w[~ 2 + Fiw\~ l + F 3 w\ + F> j +1 + F>‘ +2 
I^M]* = Hwf 2 + F>r' + F 8 w 3 + F>‘ +1 + F‘ 0 u 4 +2 
(i3s(^)J‘ = F' u w' s - 2 + F; 2 wt 1 + r ls wi + FiM +1 + F,>j +2 

ii 5 i(wi)]‘ = p>r 2 + p>r‘ + n* i + px +1 + H*> i +2 
1153 (^ 3 )]’ = piw'f 2 + p>r ! + pi*i + p<x +1 + p> 3 +2 
(LssMr = Pi>r* + p>r J + p; 3 «>j + p;x +1 + Pi>s +2 - 

and 

[tf h (Wi)) 4 = I?>r l + E\w\ + D\w * +l ' 

[//uM* = D\w'i l + D\w\ + D' 9 w' 3 +1 
(^isK)J* = D\ 2 w\~' + D[ 3 w' s + D\ 4 wi +l 
[^3i(u)i']‘ = G 2 w J" 1 + G 3 toJ + G l 4 w\ +l 
1^33(U’ 3 )]‘ = G\w\~ x + G\w\ + G' 9 w' 3 +i > 

(#35(^5)]* = G'u^r 1 + G\ 3 w\ + G\ 4 w 3 +1 
[HsiM* = H' 2 w\~ l + H 3 w\ + H\w \ +1 
{Hm(w 3 )}' = H'yw'f 1 + iliwi + P>’ +1 

= H\ 2 wf l + H\ 3 w 1 + . 

where 

e\ = c' 3 ;'/a 4 

e\ = (q 3 + c\r - c\ 2 x/2 ) /a 2 - 2 (c' 33 + ci 3 >) /a 4 

= C' n - ( 4 C ; 3 - C‘+ ,/a - C' 2 - 2 1/2 ) /a 2 + (c' 33 x + 4 C 33 + eft 1 ) /a 4 , 
E\ = (cj 3 + c\r - c£ l/i ) /a 2 - 2 (c*3 + C& 1 ) /a 4 
f; = c'+7a 4 


(D2) 


(D3) 


(D4) 
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.e^ = c*6 7 a 4 

Ej = (c; 6 + C$4 1 - C 25 ^ j A 2 -2 CJ 6 1 ^ j A 4 
E\ = C| 4 - (2C|e + 2 C 24 - Cjj 1/2 - C*- 1/2 ) /a 2 + ( C*+> + 4 C3e + C*e 1 ^ /a 4 - 
J ^9 = (^16 + ^34 1 - ^"25 ^ 2 ) /a 2 - 2 (cjg + Cjj 1 ) j A 4 
£io = c&VA 4 

£|i=C£7A 4 

E',2 = (c; 9 + cj 7 1 - C 2g ^ ^ /a 2 — 2 (cj 9 + Cj 9 j A 4 

^3 = ^7 - (2C| 9 + 2Q 7 - c‘- 1/2 - c' 2 t l/2 ) /a 2 + (cjj 1 + 4Cj 9 + cj; 1 ) /a 4 » 

E\* = (c\ g + C‘+‘- C 2g 1/2 ) /a 2 - 2 (cj 9 + Cg 1 ) /a 4 

f; 5 = c*+ 7 a 4 

FJ=C*e7A 4 

FI = (CU + C*s 1 - C£ 1/2 ) /A 2 - 2 (CU + Cje 1 ') /a 4 

^3 = ^14 — (2C34 ~' r 2^16 — ^25 ^ “ ^25 ^ ) /A 2 + (C3J 1 + 4(7^ + Cjg 1 j j A 4 • 

n = (^34 + C\V - c\t m ) /A 2 - 2 (cj 6 + a+ l ) j A a 
*s = c*+7a 4 

f‘ = c*, 7 a 4 

Fi =r (c* 6 + c- 1 - c’- 1/? ) /a 2 - 2 (c*g + c*; 1 ) /a 4 
^8 = ^44 “ (^46 - <^55 ^ ~ C55 /A 2 + (^66 * + 40^ + Cgg j A 4 > 

n = (C\e + C\t l ~ C\V' 2 ) /A 2 - 2 (c*g + C& 1 ) /A 4 
F * 0 = C'+'/A 4 

F\\ = <^69 7 A 4 

f* 2 = (c* 9 + c- 1 - c’g 1/2 ) /a 2 - 2 (c* 9 + cj 9 >) /a 4 

f; 3 = q 7 - ( 2 C‘ P + 2 Q 7 - c‘g 1/2 - c£ 1/2 ) /a 2 + (q+‘ + 4C* 9 + c*; 1 ) /a 4 > 

F\i = (^49 + ^67 1 - ^58 1/2 ) /A 2 - 2 ^Cg 9 + Cgg 1 j j A 4 

FU = cit 7 a 4 

p* = q 9 7A 4 

f* = (q 7 + c;; 1 - c’g ,/2 ) /a 2 - 2 (cj 9 + c £ l ) /a 4 

P3 = q 7 — (2Q 7 + 2CJg — C 2 g 1/2 — C 28 1 ^ 2 ^ j A 2 + ^C g9 1 + 4C 29 + C g9 j A 4 ■ 
^4 - (^37 + ^19 1 “ ^28 ^ 2 ) /a 2 - 2 ^CJ 9 + C3J 1 j j A 4 

p 5 ' = c*J7a 4 


(D5) 


(D6) 


(D7) 


(D8) 


(D9) 


(DIO) 
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*S = C£7 A 4 V 

P 7 = ( C 67 + 6 V - <^58 1/2 ) /A 2 - 2 (cj 9 + C’g ^ /A 4 

P 8 = C 47 - ( 2G 67 + 2C’g - c' + 1/2 - C* j 1/2 ) /a 2 + (c^ 1 + 4C£g + C’g /a 4 . 

p 9 = (Q 7 + G « 1 - Cs 8 1/2 ) /a 2 - 2 (cjg + Q+ 1 ) /a 4 

^0 = ^69 ! M 4 

j 

P \ 1 = ^99 1 /^ 4 

A* 2 = (cfe + C£ » - C*- 1 / 2 ) /A 2 - 2 (C8. + <£») /A 4 

P ‘ 3 = “ ( 4C 79 ~ ^88 V2 “ <?88 1/2 ) /A 2 + (C7*+» + 4Cjg + C£*) j A 4 

P 14 = ( G 79 + ^79 ‘ ~ <^88 1/2 ) /A 2 - 2 (Cjg + G^ 1 ) /A 4 

P 15 = Cm ‘/A 4 

05 = G£ ,/ 2 /A 2 

d 3 = -G'„ + (g$+ ,/2 - g;; i/2 ) /a - (g£ 1/2 + g’; 1 / 2 ) /A 2 

D* = G£ 1/2 /A 2 

^ = G^ 2 /A 2 -(G^-G* 15 )/2A' 

D8 = -G*u-(Gir /2 + G^ 2 )/A 2 • 

^ = G;r /2 /A 2 + (G-'-G* 15 )/2A 

^2 = G’; 1/2 /A 2 - (G‘7 1 - G* 18 ) j 2A j 

D; 3 = -G* 17 -(G*r /2 + G- 1 / 2 )/A 2 
^4=G*r /2 /A 2 + (G-'-G‘ 8 )/2Aj 

G 2 = G 25 1/2 /A 2 - (g‘,5 1 - GJ 4 ) j 2A | 

G 8 = -G* u - (Gjj 1/2 + Gj7 1/2 ) j A 2 1 

Gi = G’+ I/2 /A 2 + (Git , -G* 4 )/2Aj 

G« =G‘^ 2 /A 2 

Oi = ~Gi4 + (G\V /2 ~ G\V /3 ) /A - (G 8 r /2 + G 88 1/2 ) /a 2 

gS-g£ ,/8 /a 3 

G i 2 = G‘ 8 1/2 /A 2 -(Gj 7 ‘-G* 8 )/2A 

g *i 3 = -g* 7 - (g£ ,/2 + G*, 1 / 2 ) /A 2 ■ 

G 14 = G;+ 1/2 /A 2 + (g^ 1 - Gig) j 2A 


(Dll) 


(D12) 


(D13) 


(D14) 


(D15) 


(D16) 


(D17) 


(D18) 
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H\ = G‘; 1/2 /A 2 - (g\~ 6 1 - G'„) /2A | 

hj = -g\ 7 -(g‘ 2 : i/2 +g‘i 1/2 )/a 2 l 
= g£ ,/2 /A 2 + (g*+> - G* 7 ) /2A J 

^ = G*s 1/2 /A 2 - (Gig 1 - G* s7 ) /2A ' 
U's = -G* 7 - (g£ 1/2 + G* 8 1/2 ) /A 2 • 
H l g = G*+ 1/2 /A 2 + (G‘1; 1 - G’„) /2A 


(D19) 


(D20) 


*i 2 = G‘; 1/2 /A 2 

^‘13 = -C?’ 77 + (G*r /2 ~ G^ 172 ) /A - (G*+ 1/2 + G*; 1/2 ) /A 2 • 

^4 = G‘V /2 /A 2 


(D21) 


The C M and G M terms appearing in equations (D4) through (D21) are defined in appendix C by equations (C7) 
and (C32). 

Applying the difference expressions for the differential operators given in equations (D2) and (D3) to the governing 
differential equations at each finite difference station on the interior of the one-dimensional domain results in a set 
of homogeneous algebraic equations given by 


[K]{u;} = A[G]{u;} (D22) 

The matrices [ K ] and [G] are of dimension 3 (Af - 2) by 3(M + 2), where M is the number of finite difference stations. 
To reduce the system of equations further and render them symmetric, the boundary conditions are applied. The 
boundary conditions of interest in this study are given in appendix C by equations (C41) and (C42). The finite 
difference expressions for the boundary conditions are obtained directly from the central difference expressions and 
are summarized as follows: 

Simply supported edges at x = ±c 


Clamped edges at x = ±c 


wj = = 0 

w° = -wi 



) 

wl=w$* =0 

U>3 = ~ W 3 

wf +1 


> (D23) 

w\ = w? = 0 

W° = ~ W l 



1 

tu{ = wf* = 0 

ti>i = I# 2 

«,W +l : 



wj = w = 0 

«>3 = «>! 

Wj #+l ! 


(D24) 

w\ = = 0 

W? = W? 





In the boundary conditions given by equations (D24), the first derivative is defined at full stations instead of at half 
stations to simplify implementation into the computer program. 
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V 




& 


After applying the boundary conditions, the resulting system of equations constitutes a real, symmetric 
generalized eigenvalue problem given by 


*n 

*13 

*15 

( Wi ' 


911 

013 

1 

015 

*31 

*33 

} #35 

l ™3 

> = A 

931 

| 933 

} 035 

*51 

j *53 | 

s 

K W*y J 


051 

1 1 
j 953 

1 

1 955 


w x 

w z 

™5 


where 


w x = {w\ w\ 

= {«3 w l 

»T r o i 

w£ = {w£ 


1 ^ 

w?- 1 ' 
w“- l} 


Because of the self-adjoint nature of the differential operators, it follows that 


Furthermore, as a result of the bending energy being positive definite, the matrix 


#11 j #13 j #15 


! #33 ! #35 


Symmetric 


#, 


55 


(D25) 


(D26) 


- * 

II 

S 

031 = 013 1 




*5i = *r 5 

T 1 

051 = 9is | 

► 

(D27) 

’ 

*53 = *35 

053 = 035 J 




(D28) 


is positive definite. The submatrices # tJ (after assembly of the difference equations) are described in the following 
matrix equations. In each of these matrices, the ± symbol refers to the boundary condition on the loaded edges: 
the plus sign corresponds to the loaded edges being clamped, and the minus sign corresponds to the loaded edges 
being simply supported. 


Kn = 



E i 

E i 

0 

0 

0 

0 •• 

0 

0 0 

0 

0 


E $ 

E* 

E l 

0 

0 

0 •• 

0 

0 0 

0 

0 



E i 

E i 

E i 

0 

0 .. 

0 

0 0 

0 

0 




E* 

El 

El 

0 • 

0 

0 0 

0 

0 








P M- 5 
*3 

e M- 5 e M- 5 

0 

0 










E«-< 

0 









rW~3 

“3 

e M-3 

E *- 3 










pM"2 

*3 

e M-7 

. Symmetric 










(eM-Ue**- 1 ). 


(D29) 
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\( E l ±E l) 

E l 

*?0 

0 

0 

0 

0 • 

0 

0 

0 

0 

0 

0 

0 

\ 


E 7 

E l 

n 

*?o 

0 

0 

0 ■ 

0 

0 

0 

0 

0 

0 

0 




*7 

E i 

E 4 

*10 

0 

0 ■ 

0 

0 

0 

0 

0 

0 

0 



0 

E l 

E i 

E 8 

E l 

E \o 

0 • 
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0 

0 

0 

0 

0 

' 

Kl3 = 

0 

0 

0 

0 

0 

0 

0 •• 

pM-5 

*6 

e M-S 

pM - 5 

*s 

•pM- 5 
E 9 

rW-5 

*10 

0 

0 

3 

V t 


0 

0 

0 

0 

0 

0 

0 •• 

0 

pM- 4 
*6 

rM-4 

*7 

rM-4 

Eg 

pM~ 4 
E 9 

pM- 4 
E lO 

0 



0 

0 

0 

0 

0 

0 

0 • 

0 

0 

pM - 3 
*6 

E M-3 

pM—3 

*8 

n 

1 

“5 C3J 
C*3 

pM — 3 
E 10 

1 


0 

0 

0 

0 

0 

0 

0 •• 

0 

0 

0 

pM - 2 
*6 

e 7 m - j 

pM ~ 2 
*8 

E s' 2 



0 

0 

0 

0 

0 

0 

0 •• 

0 

0 

0 

0 

f M — 1 
*6 

pM— 1 
E 1 

K' l±E \ 




r(*? 3 ±*?i) 

®?4 

E h 

0 

0 

0 

0 ••• 

0 

0 

0 

0 

0 

0 

0 

- ; 


*?2 

*13 

E !< 

*?s 
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0 

0 

0 

0 

0 

0 

0 

0 

0 

s 


*?1 

*?2 

«?3 

*?4 

E U 

0 

0 

0 

0 

0 

0 

0 

0 

0 

: 


0 

*fi 

£? 2 

®?3 

E U 

*?5 

0 • • 

0 

0 

0 

0 

0 

0 

0 


K,s = 

0 

0 

0 

0 

0 

0 

0 

pM-5 

*n 

pM- 5 
*12 

pM-5 

*13 

pM- 5 
*14 

r M~ 5 
*15 

0 

0 



0 

0 

0 

0 

0 

0 

0 ••• 

0 

*11 

pM—4 

*12 

pM-4 

*13 

pM-4 

*14 

•rM -4 
*15 

0 



0 

0 

0 

0 

0 

0 

0 

0 

0 

pM-3 

*11 

pM—3 

*12 

pM-3 

*13 

pM-3 

*14 

pM-3 

*15 



0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

pW — 2 
*11 

pM-2 

*12 

pM — 2 
*13 

r Af-2 

*14 



0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

ip Af — 1 
*11 

jrM-1 

*12 



(f|±f 6 3) f* f» 0 0 0 0 0 

F| F 9 3 F 3 , 0 0 0 

F* F* Fj 0 0 0 

F| F| F* 0 0 


K 33 = 


* Symmetric 


?A/"5 pM-5 


7 »A/ — 4 nA/-4 


?A/— 3 pA/“3 


pM -3 

*10 

jj .M -2 

*Q 
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K35 = 




p2 

M5 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

F h 

*?3 

F 3 
” 1 4 

*?s 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

F U 

F n 

p 4 
M3 

F U 

rfs 

0 

0 

u 
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0 

0 

0 

0 

0 

0 

A' 5 

Ml 

pb 
1 12 

F h 

r?4 

*?5 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

nM-5 

" Ml 

pM-5 

M2 

pM-5 

M3 

pW-5 

M4 

pM-b 

M5 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

rM- 4 
Ml 

pM-4 
* 12 

pM-4 

M3 

pM-4 

M4 

pM-4 

M5 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

pM-3 

Ml 

pM- 3 
* 12 

pM- 3 
r 13 

pM-3 

M4 

pM-3 

M5 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

pM — 2 
Ml 

pM~2 

M2 

pM-2 

M3 

pM-2 
* 14 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

pM - 1 

Ml 

pM — 1 

M2 
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K55 = 


jiPfl) P?4 P? 5 


p3 p3 p3 
M3 M4 M5 

p4 p4 p4 
M3 M4 r 15 

^3 P 14 


0 

0 


0 
0 

0 0 


p5 n 
14 M5 U 


0 

0 

0 

0 

pM- 5 
13 


L Symmetric 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

pAf— 5 
M4 

pAf-5 

M5 

0 

0 

pAf— 4 
M3 

pAI— 4 
M4 

pM-4 

M5 

0 


pM-3 

M3 

p! f — 3 
M4 

pA/— 3 
M5 



nM-2 

M3 

pM-2 


(D34) 


The elements of the matrices are defined by equations (D4) through (D12) and M is the number of difference 
stations. The gij submatrices are similar to the K%j submatrices with the Ej being replaced by the F? being 
replaced by Gf, and the Pf being replaced by the H? (see eqs. (D13) through (D21)). 
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Appendix E 

Modifications for the Stress Loading Case 

When specifying the stress loading case, differential equations (B57) and boundary conditions (B62) must be 
modified to eliminate any rigid body displacements of the plate. The modifications used in this study are presented 
in this appendix. 

Examination of equations (B57) indicates that the first and eighth differential equations are independent of the 
others. From the eighth differential equation, t/g = 0, the solution can be obtained directly as y$ = t/ 88 , a constant. 
To separate the system and make use of PASVAR, the equations are redefined as follows: 


5 

II 

§ 




3/5 = ye 

2/9 = yn 

yi = yi 




£ 

II 

s 

yio = yi2 

<ci 

M 

II 

w 




£7 = y9 

2/11 = yi3 

2/3 = y4 



2/4 = Vb 

ys = yio 

2/12 = yi4 


(El) 


Redefining the unknown variables leads to the following set of differential equations for PASVAR: 


^12 £23/88 + ?52/7 + $62/8 + £72/9 

Vi = --7—1V4 + -r—f . — : 

An AnDenl 

~l 3i4j2 - ?32/88 + ?6j/7 + ?82/8 + $92/9 

2/2 = — A — 72/5 + 


2/3 


y '4 = 72/1 + 
V 5 - 372/2 + 
2/6 = 372/3 + 


^11 

5j4i2 _ 
-72/6 + 


.AnDenl 

?42/88 + f72/7 + ?92/8 + fl02/9 


An "" ' i4nDenl 

0152/10 + 0162/11 + 0172/12 


A 6 6Den2 

0162/10 + 0182/11 + 0192/12 
,4 66 Den2 

0172/10 + 0192/11 + 0202/12 


A6eDen2 


(E2a) 

(E2b) 

(E2c) 

(E2d) 

(E2e) 

(E2f) 
(E2g) 
(E2h) 
(E2i) 

(E2j) 

(E2k) 

i\-2U“u ' 1 v v /•mowo; ■ A iw (E21) 

11 An 

Upon solution of this system of equations, the displacement component yo is obtained by direct integration of the 
following differential equation: 


2/7 = “72/io 
Vs = -37$! 1 

y'g =• -57S12 
_ ^11-^22 

VlO “ 


‘11 


4 7O52/4 + 37n 4 y 5 + Jh^nye) + -7^72/7 

An 


yii = 


Via ~ 


— — — Y — 37(03^0 4* 7^42/4 4* 37002/5 + 57^112/6) 4- 

An 

AnAaa - A? 


3An _ 
1 — iy* 

An 
5Ai 


“ 57 (^ 10^0 + 7^122/4 + 370112/5 4* 57 O 132 / 6 ) + “7^72/9 


^ , ?iy88 + f2y7 + ?3ys + j 4 y9 

y 0 = — vo 4 * t — ~ — 

An An Deni 


(E3) 


35 


APPENDIX E 


The expressions for the stress resultants given by equations (B63) become 


\f yes + tii/7 + {3V8 + UV9 

10 Deni 

(E4a) 

ivr _ ^22/88 + &V7 + tei/8 + {7$9 

xl Deni 

(E4b) 

nr $ZV88 + tei/7 + f8 j'S + fyyg 

N ’ a ~ Deni 

(E4c) 

XT <42/88 + <72/7 + <92/8 + <102/9 

N “ = Deni 

(E4d) 

n„. a » a »- a » V0 + ±i N x , 0 

An An 

(E4e) 

» r M 1 M 2 — A}2 ~ , A 12 hr 

Myl= 1V4+ A Nsl 

An >ni 

(E4f) 

\r 0 ^ 11^22 “^12 - , -^12 kt 

Ny3=Z . 73/5+. N x3 

a 11 All 

(E4g) 

AT r ^11^22 ^12 *- . A 12 .. 

NyS — 5 , 7V6 + . N xi 

Ml A 11 

(E4h) 

w OlbVlQ + 0162/11 + 0172/12 

"**' “ Den2 

(E4i) 

AT 0162/10 + 0182/11 + 0192/12 

= Den2 

(E4j) 

xr 0173/10 + 019^11 4- 0202/12 

xyi ~ Den2 

(E4k) 

The boundary conditions for the system of differential equations (E2) are obtained using plate equilibrium 
This procedure is described as follows. Using the expressions for the stress resultants given by equations (B39) 
through (B41) with equations (B43) and (El) gives 

ye 8 — N x0 n t 4* Jv*in 2 + + N x sttio 

(E5) 

As mentioned previously, y&8 is a constant, and this expression represents the axial force divided by 2A at any croos 
section of the plate. Thus, equating the applied force at the loaded edges with equation (E5) multiplied by 2A yields 

V88 = ~bN Q x 

(E6) 


The 14 boundary conditions corresponding to the uniform compressive stress loading, given by equations (B32), are 


N x0 (±c) = -N° 
W*i(±c) = 0 
N x3 (±c) = 0 
W*s(±c) = 0 


Nxyl (it) — 0 I 
N xy3 (±c) = 0 1 
JV* v5 (±c) = 0 


(E7) 


Tho corresponding 12 boundary conditions for differential equations (E2) are obtained by eliminating two boundary 
conditions from equations (E7). The boundary conditions to be eliminated are obtained directly from equation (E5) 
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by specifying 


N x0 (±c) = 


Hi(±c) 


JV sl (±e) = 0 
JV« 3 (±c) = 0 

Equilibrium equation (E5) requires that N x s(±c) = 0, since Qi 0 (±c) is nonzero. Thus, the boundary conditions of 
differential equations (E2) are 

02/88 + 02/7 + 6*2/8 + 02/9 = — iV^Denl 
$ 22/88 + $ 52/7 + $ 62/8 + 02/9 = 0 


$ 32/88 + 02/7 + 02/8 + 02/9 — 0 

► {x = ±c) (E9) 

0153/10 + 0162/11 + 0172/12 = 0 
0162/10 + 0182/11 + 0191/12 == 0 
0172/10 + 0192/11 + 0202/12 = 0 

The remaining step in solving the stress loading case is to determine the value of the constant Vo- The correct 
value of vo is defined to be the value which makes the residual normal force acting on the unloaded edges negligible 
compared with the applied loading. The net residual normal force on the unloaded edges is 


The applied load is 


r+c 

P y (±b) = A J \Ai 2 u' 0 (x) 4- A 22 V 0 ] dx 
P x (±c) = 2 f N x (± c, y) dy = 2Ai/ 88 (±c) 


The constant v 0 is related to the generalized displacements by 

«o - ~^^[«°(c) - Uo(-c)) (E12) 

When the pl»*e is .tress loaded, the displacements of the loaded edges, u°(±c), are not known a priori. However, 
the prebuckling problem is linear and vq can be obtained by interpolation. 
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Appendix F 

User’s Guide and Sample Problems 

For a constant cutout size and plate width, the pro- 
gram can calculate the buckling load and mode shape 
for a plate for several aspect ratios. The buckling results 
for the plate having the initial half-length, STARTC, 
are calculated. Then, the plate half-length is increased 
by CINC, and the buckling results for the plate with the 
new aspect ratio are calculated. This incremental pro- 
cess continues until the buckling results for the plate 
having the final half-length, STOPC, are calculated. 
The program input data are simple and are contained 
on four card images. The data on each card image are 
described as follows: 

Card Is CONV, CONF, NSTATS 

CONV, CONF convergence tolerances on the pre- 
buckling displacements and stresses, 
respectively 

NSTATb number of finite difference stations 

The values cf CONV and CONF represent the final 
maximum errors in the displacements and stresses cal- 
culated by PASVAF. These values do not have to be 
know n exactly and should be specified in accordance 
with the magnitudes of the displacements and stresses 
expected. CONV is much smaller than CONF since the 
displacements are several orders of magnitude smaller 
than the stresses. 

Card 2: A, XI, B, STARTC, CINC, STOPC, IPRINT, 
IBCON, ICTOUT, ILOAD 

A, XI cutout dimensions (see fig. 1) 

B plate half-width (see b on fig. 1) 

STARTC initial value of the plate half-length 
(see c on fig. 1) 

CINC plate half-length increment for 

multiple aspect ratio calculations 

STOPC terminating value of the plate half- 

length 

IPRINT output option: 

1 minimum program output 

2 maximum program output 


IBCON boundary condition option: 

1 loaded edges clamped 

2 loaded edges simply supported 

ICTOU T cutout shape option. 

1 rectangular cutout given by 
x — A and y = XI 

2 elliptical cutout given by 
(x/A)**2 + (y/XI)**2 — 1 

ILOAD loading option: 

1 uniform edge displacement 

2 uniform edge stress 

Card 3: All, A12, A22, A66 

A11,A12, membrane stiffness coefficients for 
A22,A66 orthotropic plates 

Card 4: Dll, D12, D22, D66 

D11,D12, bending 3tiffness coefficients for 
D22.D66 orthotropic plates 

To execute multiple analyses, repeat cards 2, 3, and 
4 for each additional problem. The program execution 
halts on reaching the end of input data. 

To illustrate program use, two sample problems are 
presented. 

Sample ProbL n I 

The buckling load, buckling coefficient, and critical 
end shortening for a square orthotropic plate with a 
centrally located circular cutout are calculated. The 
plate is 10 inches on a side, and multiple analyses 
are executed for ratios of cutout diameter to plate 
width from 0 to 0.6. The stiffness coefficients of the 
orthotropic plate are given by 

‘18.698 0.566 0 

[A tj ) = 0.566 1.617 0 x 10 5 Ib/in. 

[ 0 0 0.832 

‘ 15.582 0.472 0 ] 

[Dij] = 0.472 1.348 0 x 10 2 lb/in. 

[ 0 0 0.693 

All edges of the plates are simply supported. The 
loading is a uniform edge displacement. The corre- 
sponding input data are 
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APPENDIX F 


Card 1: 1.0 6.2E+5 62 
Card 2: 0.0 0.0 5.0 5.0 0.0 5.0 1 2 2 1 
Card 3: 18.698E+5 0.566E+5 1.617E+5 0.832E+5 
Card 4: 15.582F+2 0.472E+2 1.348E+2 0.693E+2 

Card 5: 0.5 0.5 5.0 5.0 0.0 5.0 1 2 2 1 
Card 6: 18.698E+5 0.566E+5 1.617E+6 0.832E+5 
Card 7: 15.682E+2 0.472E+2 1.348E+2 0.693E+2 

Card 8: 1.0 1.0 5.0 5.0 0.0 5.01221 
C?'d 9: 18 . 698E+5 0.566E+5 1.617E+5 0.832E+6 
Card 10: 15.582E+2 0.472E+2 1.348E+2 0.693E+2 


Card 11: 1.5 1.5 5.0 5.0 0.0 5.0 1 2 2 1 
Card 12: 18.698E+5 0.566E+5 1.617E+5 0.832E+5 
Card 13: 15.582E+2 0.472E+2 1.348E+2 0.693E+2 

Card 14: 2.0 2.0 6.0 5.0 0.0 5.0 1 2 2 1 
Card 15: 18.698E+5 0.566E+5 1.617E+5 0.832E+5 
Card 16: 15.582E+2 0.472E+2 1.348E+2 0.693E+2 


Card 17: 2.5 2.5 5.0 5.0 0.0 5.0 1 2 2 1 
Card 18: 18.698E+5 0.566E+5 1.617E+5 0.832E+5 
Card 19: 15.582E+2 0.472E+2 1.348E+2 0.693E+2 


Card 20: 3.0 3.0 6.0 5.0 0.0 5.0 1 2 2 1 
Card 21: 18.698E+5 0.566E+5 1.617E+5 0 832E+5 
Card 22: 15.582F+2 0.472E+2 1.348E+2 0.693E+2 


As described above, cards 2, 3, and 4 are repeated 
for each cutout size. The results from BUCKO are 
summarized for ratios of the hole diameter to plate 
width d/W below. 


d/W 

Buckling 
load, lb 

Buckling 

coefficient 

Critical end 
shortening, in. 

0 

2037 

4.50 

0.00110 

.1 

1967 

4.37 

.00108 

.2 

1777 

3.93 

.00106 

.3 

1641 

3.63 

.00114 

.4 

1560 

3.45 

.00126 

.5 

1512 

3.34 

.00145 

.6 

1564 

3.46 

.00186 


Sample Problem 2 

The buckling load, buckling coefficient, and criti- 
cal end shortening for a rectangular orthotropic plate 
with a centrally located rectangular cutout are calcu- 
lated. Buckling results are calculated for plates having 
a width of 10 inches and lengths of 20 inches, 25 inches, 
and 30 inches. The rectangular cutout is 2 inches wide 
and 4 inches long with the 2-inch-wide side perpendic- 
ular to the loading. The stiffness coefficients used in 
sample problem 1 are also used for this sample prob- 
lem. The loaded edges are clamped, the unloaded edges 
are simply supported, and the loading is a uniform edge 
stress. The corresponding input data arc 


Card 1: 1.0 6.2E+5 42 

Card 2: 2.0 1.0 5.0 10.0 2.5 15.0 1112 
Card 3: 18.698E+5 0.566E+5 1.617E+5 0.832E+5 
Card 4: 15.582E+2 0.472E+2 1.348E+2 0.693E+2 


The results obtained from BUCKO are summarized 
as follows: 


Aspect 

ratio 

Buckling 
load, lb 

Buckling 

coefficient 

Critical end 
shortening, in. 

2.0 

1997 

4.37 

0.00283 

2.5 

1729 

3.82 

.00280 

3.0 

1511 

3.34 

.00283 
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